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ABSTRACT 


Recent  advances  in  thermal  imaging  technology  have  resulted  in  the  fielding  of 
two-dimensional  array  detector  based  imaging  systems.  These  designs  have  been  labeled 
second-generation,  and  are  rapidly  replacing  first  generation  systems  having  linear 
detector  arrays  with  a  parallel  scan  type  architecture.  It  has  been  postulated  that  first 
generation  prediction  models  are  not  applicable  to  second  generation  systems.  In 
particular,  the  minimum  resolvable  temperature  difference  (MRTD)  modeling  needs 
refinement  in  the  areas  of  sampling,  quantization  noise,  and  array  non-uniformities  in 
order  for  it  to  be  applied  to  second  generation  systems.  The  present  industry  standard  for 
MRTD  is  the  Night  Vision  FLIR92  Model.  Results  from  the  FLIR92  Model  and  the  two 
well  known  first  generation  models  will  be  presented  and  compared  with  experimental 
measurements  made  on  two  thermal  imaging  systems  available  at  the  Naval  Postgraduate 
School. 
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I.  INTRODUCTION 


A.  REVIEW  OF  LITERATURE 

Advances  in  thermal  imaging  technology  have  resulted  in  the  fielding  of 
two-dimensional  array  detector  based  imaging  systems.  This  advancement  in  technology 
resulted  in  a  reexamination  of  the  applicability  of  the  U.S.  Army  Night  Vision 
Laboratory's  1975  Static  Performance  Model.  This  model  was  updated  as  FLIR90  in  1990 
and  finally,  it  was  updated  to  FLIR92  to  predict  the  performance  of  both  first  and  second 
generation  TISs  [Ref.  9]. 

The  minimum  resolvable  temperature  difference  (MRTD)  is  generally  considered 
to  be  a  very  good  performance  indicator  for  TISs.  The  MRTD  is  experimentally 
measurable  and  it  has  been  subject  for  predictive  modeling.  A  good  background  reference 
to  appreciate  the  significant  issues  in  MRTD  modeling  is  provided  by  Lloyd’s  classic  text 
on  thermal  imaging  systems  [Ref.  1].  The  MRTD  model  proposed  in  Lloyd’s  book  is  a 
simplified  treatment  of  the  earlier  work  of  Ratches  [  Ref. 4]  based  on  the  eye-brain  match 
filter  model. The  Ratches  MRTD  model  also  known  as  the  Ratches-Lawson  model  was 
based  on  the  matched  filter  and  synchronous  integrator  model  [Ref.  6].  The  Vortman  and 
Bar-Lev  model  [ref.  5]  has  also  been  proposed  to  improve  the  MRTD  predictions  and 
compared  to  the  MRTD  models.  However,  questions  still  remain  concerning  the  accuracy 
of  these  models  [Ref.l],  [Ref.  4],  [Ref.  5],  [Ref.  6].  The  most  important  issues  in  MRTD 
modeling  are; 

1.  Visual  Perception 

Much  research  has  been  done  on  this  area  and  different  theories  and  models  have 
been  developed.  Most  of  this  work  tries  to  model  the  MTF  of  the  eye-brain  combination 
[Ref.  4,  Ref.  8].  Recently  a  simple  eye  MTF  model  which  assumes  “a  linear  dependence 
up  to  a  normalized  peak  response,  at  approximately  0.4  cyc./mr. ,  followed  by  a  decay  in 
response  predicted  by  the  optical  in-focus  OTF  associated  with  the  finite  pupil  size  of  the 
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eye.”  has  been  applied  to  MRTD  modeling  by  Pieper  [Ref.  8.  P.  263].  Figure  I.l  shows 
this  bandpass  MTF  model  for  eye  filter.  A  theory  for  eye  and  brain  processing  to  detect 
sinusoidal  gratings  with  fluctuations,  which  are  caused  by  video  noise,  was  presented  by 
Schnitzler  [Ref  33].  The  eye-brain  combination  is  modeled  as  two  filters  by  Komfeld 
and  Lawson  [Ref. 34].  The  spatial-frequency  filtering  by  the  eye  near  to  the  sine-wave 
modulation  threshold  is  discussed  by  Schnitzler  [Ref. 35].  The  eye  and  brain  are 
considered  as  two  different  filters  and  compared  by  Campbell  [Ref  36]. 


2.  Sampling 

Image  sampling  is  the  processing  between  the  continuous  scene  (input  image), 
discrete  sampled  image  signal  and  continuous  reconstructed  image.  Signal  conversions 
between  them  caiase  degradation,  aliasing  and  bliirring,  in  TISs.  The  NVEOD  FLIR92 
TISs  Performance  Model  and  its  noise  sampling  effects  is  introduced  by  D'  Agostino  and 
Scott  [Ref.  16].  The  calculation  of  magnitude  of  the  image  degradation,  related  with 
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aliasing  and  blurring  is  explained  by  Huck,  Halyo,  and  Park  [Ref.  17].  Some  methods  to 
reduce  the  effects  of  the  aliasing  are  proposed  and  the  application  of  naicroscanning  is 
explained  in  Watson,  Muse,  and  Blommel,  [Ref.  18].  A  method  to  increase  the  sampling 
rate  for  undersampled  staring  systems  is  proposed  by  Chambliss,  Dawson,  and  Borg 
[Ref.  19].  The  effects  of  the  sensor  aliasing  to  the  target  detectability  is  showed  by 
Meitzler  and  Gerhart  [Ref.  30].  The  MTF  of  the  system  due  to  effects  of  the  image 
processing  is  analyzed  by  Park,  Schowengerdt,  and  Kaczyniski  [Ref.  31]. 

3.  Noise  Characteristics 

Detector  noise  is  the  dominant  noise  for  the  first  generation  TISs.  Improvements 
in  the  TISs  introduced  the  3-D  noise  concept,  because  modem  systems’  noise  show  more 
directionality  than  the  first  generation  TISs  and  this  concept  was  used  for  FLIR92.  3-D 
noise  methodology  and  its  noise  components  are  introduced  by  D'Agostino.  3-D  noise 
and  the  usage  of  directional  averaging  operators  (D  operators)  are  defined  by  Webb  and 
Bell  [Ref.  12].  The  3-D  noise  elements  and  the  noise  correction  factors  are  explained  by 
Scott  and  D'Agostino  [Ref  13].  Each  of  the  3-D  noise  components  is  presented  by  Webb 
[Ref  14]. 

Laboratory  measurements  have  been  done  objectively  to  compare  the  TISs’ 
measured  MRTD  and  predicted  MRTD  by  models.  The  usage  of  automatic  target 
recognition  units  to  eliminate  the  human  observer's  subjectivity  is  being  tested  for 
different  directions.  Williams  [Ref.  21]  describes  a  technique  to  replace  the  observer  by  a 
CCD  array  camera  and  computer.  Williams  [Ref.  22]  explains  the  requirements  of 
objective  measurement.  The  measurement  techniques  and  requirements  to  standardize  the 
MRTD  results  between  the  laboratories  have  also  been  given  attention  [Ref  23].  Recently 
a  report  has  been  disseminated  which  covers  the  proceedings  of  the  1990  NVEOD  FLIR 
modeling  workshop  [Ref.  29]. 
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B.  OVERVIEW  OF  THE  THESIS 


This  chapter  reviews  the  fundamentals  and  basic  parameters  of  thermal  imaging 
systems  and  briefly  mentions  the  Johnson  criterion. 

Chapter  n  defines  the  Batches  Model  and  the  Lloyd  Model  which  are  applicable 
to  first  generation  non-staring  systems.  It  can  be  shown  that  these  two  models  can  be 
related  and  are  equivalent  to  each  other  at  low  frequencies.  Also,  in  this  chapter  the  U.S. 
Army  Night  Vision  Laboratory's  FLIR92  TIS  Performance  Model  for  MRTD  is  compared 
with  the  Batches  and  the  Lloyd  models. 

Chapter  IE  presents  the  FLIR92  TIS  Performance  Model  under  three  parts.  In  the 
first  part,  it  briefly  shows  the  MTFs,  which  are  used  by  this  model.  In  part  two,  it 
describes  the  3-D  noise  concept,  noise  applications  and  noise-data  input  groups.  In  part 
three,  it  defines  the  Minimum  Resolvable  Temperature  Difference  (MRTD)  and 
Minimum  Detectable  Temperature  Difference  (MDTD)  by  using  3-D  noise  concept. 

Chapter  IV  covers  the  effect  of  aliasing,  which  is  not  covered  in  FT  .IRQ?  It 
mathematically  models  the  aliasing  effect  including  phase  artifacts  due  to  relative  shifts 
between  the  four-bar  image  and  detector  array.  Computer  simulations  of  results  are 
shown  to  clarify  the  significance  of  this  effect.  The  microscanning  for  reducing  the 
aliasing  effect  is  discussed  [Ref.  18]. 

In  Chapter  V  the  laboratory  measurements,  which  were  done  with  Amber 
engineering  AE4128  IR  Imaging  System  and  Mitsubishi  Electronics  IR-M500  Imager  are 
presented.  Both  experimental  MRTD  measurements  and  objective  measurements  were 
done  during  the  laboratoiy  work. 

In  Chapter  VI  the  comparisons  of  the  FLIR92  MRTD  to  the  measured  MRTDs, 
which  were  presented  in  Chapter  V,  are  presented.  Also,  in  this  chapter  the  limitations  of 
the  FLIR92  model  are  presented. 

Appendix  A  gives  the  parameters  of  the  systems,  which  are  used  in  this  thesis. 
Appendix  B  provides  the  MTF  equations  of  FLIR92,  which  are  listed  in  Chapter  HI. 


4 


Appendix  C  presents  the  data  for  experimental  MRTD  and  objective  measurements. 
Appendix  D  is  short-listing  outputs  of  the  FLIR92  for  Amber  engineering  AE4128  IR 
Imaging  System  and  Mitsubishi  Electronics  IR-M5(X)  Imager. 

A  table  of  symbols,  which  shows  the  symbols,  acronyms  and  abbreviations,  is 
included  at  the  beginning  of  the  thesis. 

C.  FUNDAMENTALS  OF  THERMAL  IMAGING  SYSTEMS 

A  Thermal  Imaging  System  collects  the  radiated  infrared  energy  of  target  and 
background,  and  changes  into  a  visible  image  for  an  observer  to  detect,  recognize,  and 
identify  the  target.  There  must  be  a  large  enough  temperature  difference  between  the 
background  and  target  to  overcome  the  noise  sources,  nonuniform  background 
temperature,  and  atmospheric  attenuation.  With  both  scanning  mirrors  Figure  1.2 
represents  a  serial  scanning  TIS  [Ref.  1]. 


DETECTOR 

PROJECTION 


Figure  1.2.  A  Simplified  Scanning  TIS  ( From  Ref.  1 ) 
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However,  with  only  one  scanner  the  same  system  would  represent  the  common 
module  FLIR  with  parallel  scanning  architecture  [Ref.  1].  Lastly,  with  direct  imaging  on 
to  the  two  dimensional  detector  array,  i.e.  no  mechanical  scanning,  the  system  would 
represent  a  staring  second  generation  design  which  is  currently  of  interest  [Ref.  1], 

[Ref.  9]. 

D.  BASIC  PARAMETERS  OF  THERMAL  IMAGING  SYSTEMS 

TISs  have  different  parameters  that  describe  different  facets  of  the  system.  In  this 
section  we  are  going  to  review  definitions  for  Noise  Equivalent  Temperature  Difference 
(NETD),  Detector  Angular  Subtense  (DAS),  Optical  Transfer  Function  (OTF),  Minimum 
Resolvable  Temperature  Difference  (MRTD),  and  Minimum  Detectable  Temperature 
Difference  (MDTD). 

1.  NETD 


This  was  the  most  used  parameter  to  define  the  TIS's  performance  to  detect  small 
signals  in  noise.  The  following  definition  shows  the  proper  way  to  explain  this  parameter: 
"  the  NETD  is  the  blackbody  target-to-background  temperature  difference  in  a  standard 
test  pattern  which  produces  a  peak-signal  to  rms-noise  ratio  (  SNR  )  of  one  at  the  output 
of  a  reference  electronic  filter  when  the  system  views  the  test  pattern"  [Ref  3,  p.  166] 

If  we  define  the  temperature  difference  between  target  and  background  as  AT, 
signal  voltage  as  V^,  and  noise  voltage  as  V„,  can  be  written  in  the  following  form, 
by  using  the  NETD  definition  [Ref  1]: 

_  AT 

V„  NETD  (1.1) 
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From  Equation  (LI)  we  can  obtain  NETD: 


NETD  =  ^  (1.2) 

V  n 

As  can  be  seen  from  Equations  (LI)  and  (1.2)  when  the  signal  voltage  is  equal  to  the  noise 
voltage,  the  NETD  will  be  equal  to  the  temperature  difference  between  target  and 
background. 


2.  DAS 

The  detector  collects  the  radiation  in  the  limits  of  its  Field  of  View  (FOV).  We 

can  define  FOV  as  Horizontal  Field  of  View  (HFOV)  and  Vertical  Field  of  View 

(VFOV).  Figure  1.3  shows  the  DAS  and  FOV  of  a  TIS.  In  this  figure  (X  is  Horizontal 
DAS,  (3  is  Vertical  DAS. 

HFOV 


Vertical 

DAS 

VFOV 


Figure  1.3.  Detector  Angular  Substense  ( From  Ref.  1) 

By  using  small  angle  approximation,  DAS  can  be  defined  in  the  following  form 
[Ref  2,  p.82]: 

DAS=(detector  size)  /  ( effective  focal  length )  (1.3) 


FLIR 
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3.  OTF 


OTF  is  the  main  parameter  specifying  the  frequency  response  of  an  electro-optical 
system.  It  is  composed  of  two  parts  described  by  [Ref  1]. 

OTF=MTFejPTP  a  4) 


where  MTF  stands  for  magnitude  transfer  function  and  PTF  stands  for  phase  transfer 
function.  The  MTF  is  also  referred  as  the  modulation  transfer  function  which  can  be 
measured  according  to  the  following  guidelines.  By  using  Figure  1.4,  the  modulation  can 
be  explained  in  the  following  way:  "Modulation  is  the  variation  of  a  sinusoidal  signal 
about  its  average  value. "  [Ref  2,  p.65] 


Modulation  =  M  = 


B  max-B  min 
B  max+B  min 


AC 

DC 


(1.5) 


Equation  (1.5)  leads  us  to  the  calculation  of  MTF, 


MTF= 


Mi: 


(1.6) 
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Figure  1.5  shows  a  simple  block  diagram  of  a  TIS,  as  we  can  see  from  this  figure  TIS  has 
four  main  blocks  [Ref.  1].  The  system's  overall  MTF  is  the  product  of  the  MTF  values  of 
these  blocks. 


Optics  Detector  Electronics  Display 


Figure  1.5.  Simple  Block  Diagram  of  a  TIS 

MTFs  of  this  block  diagram  are  defined  as  follows: 

a.  Optical  diffraction  MTF  ( H^,): 

Hopt  — i  acos(^)  —  (— )y  1  —  (— )  (1.7) 

is  the  optics  aperture  diameter,  X  is  the  wavelength  of  light,  and  f,  is  the  spatial 
frequency  in  the  horizontal  direction. 

D=X  foe  (1.8) 

defines  an  optic  cutoff  frequency,  4^;  if  we  insert  Equation  (1.8)  into  Equation  (1.7),  we 
can  write  H„p,  in  the  following  way: 

Hopt  =  ll"  acos(^)  -  (1.9) 

loc  loc  Y  lOC 

b.  Detector  spatial  filter  ( H^^,): 


_  sin(7ifx  g) 
ndet  - 
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(1.10) 


Detector  cutoff  frequency  is, 


fdstc  —  (X  (Ill) 

By  using  Equation  (1. 10)  and  (Ill)  we  will  have  the  following  equation,  which  shows 
>  will  go  to  ^^fo  when  fx=  f^etc  • 

/f'iet- Sinc(“)  (1.12) 

< .  I'ilectronics  ( Hrfect)^ 

As  a  tunction  of  temporal  frequency,  MTF  is 

=  a  13) 

and  PTF  is, 

/>7F=arcfan(-^)  (,.,4) 

Heiect  can  be  defined  as  a  function  of  spatial  frequency  by  using  Equation  1. 15  for  parallel 
scanning  systems. 


(1.15) 


Vx  is  the  angular  scanning  velocity.  H.!,..,  as  a  function  of  spatial  frequency,  is  usually 
defined  in  this  way: 


(1.16) 
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and  PTF  is, 


PTF  =  arctan(-^ )  (1. 1 7) 

d.  Display  ( : 

Equation  1.18  shows  spot  size  limited  display  MTF,  which  is  a  gaussian 
distribution  and  b  is  the  CRT  spot  size  parameter  in  this  equation.  This  is  equivalent  to 
assuming  that  the  display  can  be  described  by  an  incoherent  point  spread  function  with  a 
gaussian  space  dependence  [Ref.  38]. 

Hdisp  ~  ^  (1-18) 

The  equivalent  of  this  b  value  is  presented  in  Equation  I.  19  [Ref.  1]. 

h  =  2n‘^<j^  (1.19) 

Typical  ^  value  is  between  0.125  and  0.5 

4.  MRTD 

MRTD  is  the  minimum  temperature  difference  between  the  background  and 
standard  four  bar  pattern  when  four  bars  can  be  resolved  by  an  observer.  MRTD  is 
generally  considered  to  be  a  very  good  performance  measurement  for  TISs.  Figure  1.6 
shows  the  standard  four  bar  pattern  for  four  different  spatial  frequencies.  In  Figure  1.6  , 
W,  the  bar  width  is 

W  =  (1.20) 
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There  is  a  required  7:1  aspect  ratio  between  the  length  and  width  of  the  bar  patterns. 
Equation  1.16  and  Figure  1.6  prove  that  spatial  frequency  increases  as  target  size  gets 
smaller,  which  explains  why  the  pattern  as  a  whole  has  to  shrink  as  fj.  increases.  A 
detailed  description  of  the  MRTD  by  Ratches  [Ref.  4]  and  Lloyd  [Ref.  1]  is  going  to  be 
mentioned,  in  the  next  chapter. 


Figure  L6.  Four  Bnr  Patterns 


5.  MDTD 


MDTD  is  the  minimum  temperature  difference  between  background  and  a  square 
target  which  can  be  detected  by  an  observer.  Figure  1.7  shows  this  reference  target 


1/f 

I - 1 


Figure  1.7.  Reference  Square  Target 

The  MDTD  is  not  applied  as  a  performance  measure  of  TISs  as  frequently  as  the 
MRTD  and  it  will  not  be  discussed  in  the  thesis. 


12 


E.  JOHNSON  CRITERION 


Johnson  used  eight  military  vehicles  and  one  observer  to  develop  his  visual 
discrimination  methodology.  During  the  experiment,  he  placed  bar  charts,  whose  length 
is  equal  to  the  target  and  width  is  the  critical  target  dimension,  between  the  observer  and 
the  targets.  In  classifying  his  laboratory  results  he  divided  the  visual  discrimination  into 
the  following  four  levels  and  listed  the  required  square  wave  spatial  frequency  resolution 
requirement  for  each. 

1 .  detection 

2.  orientation 

3.  recognition 

4.  identification 

These  discrimination  levels ,  their  meanings,  and  their  required  resolution  cycles 
across  the  minimum  dimension  are  shown  in  Table  I.l.  [  Ref.  2,  p.415  ] 

Table  LI.  Johnson’s  Results 
(  From  Ref.  2  ) 


DISCRIMINATION 

LEVEL 

MEANING 

CYCLES  ACROSS 
MINIMUM 
DIMENSION 

Detection 

An  object  is  present 
(object  versus  noise) 

1.0  ±0.025 

Orientation 

The  object  is  approximately 
symmetrical  or  unsymmetrical 
and  its  orientation  may  be 
discerned  (  side  view  versus 
front  view). 

1.4±0.35 

Recognition 

The  class  to  which  the  object 
belongs  (e.g.,  tank,  track,  man) 

4.010.80 

Identification 

The  object  is  discerned  with 
sufficient  clarity  to  specify  the 
type  (e.g.,  T-52  tank,  friendly 
jeep) 

6.411.50 

Johnson's  methodology  provides  50%  probability  of  discrimination.  By  using  the 
Johnson  criterion,  we  can  write  an  equation,  which  shows  the  connection  between  the 
target  MRTD  and  model  ( four  bar  )  MRTD  [Ref  2],  as  follows: 

MRTD,^  „  =  MRTD  (1,21) 

©  is  the  aspect  ratio  of  the  target's  size,  the  ratio  of  the  maximum  target  dimension  to  the 
minimum  target  dimension,  N  is  the  required  number  of  cycles  for  detection,  is  the 

aspect  ratio  of  the  four  bar  pattern.  In  the  case  of  interest  Taspect  =7  . 

F.  MOTIVATION  FOR  TfflS  WORK 

Recent  advances  in  thermal  imaging  technology  have  resulted  in  the  fielding  of 
two-dimensional  array  detector  based  imaging  systems.  These  "staring"  designs  have 
been  labeled  second-generation  and  are  rapidly  replacing  existing  first  generation 
systems,  some  with  single  detectors  but  primarily  one  dimensional  arrays  of  detectors 
having  parallel  scan  type  detection.  It  has  been  postu’-^ted  that  first  generation  prediction 
models  are  not  applicable  to  second  generation  systen., .  The  MRTD  modeling  process 
needed  refinement  in  the  areas  of  sampling,  quantization  noise,  and  array 
non-uniformities  in  order  for  it  to  be  applied  to  second  generation  systems. 

The  purpose  of  this  thesis  is  to  compare  the  FLIR92  Performance  Model  for  TISs 
to  the  Ratches  MRTD  Model  and  the  Lloyd  MRTD  Model  for  first  generation  TISs  and 
to  measured  data  for  second  generation  TISs. 


14 


II.  FIRST  GENERATION  THERMAL  IMAGING  SYSTEMS 


A.  SAMPLE  SYSTEM 

In  this  chapter,  the  Ratches  Model  and  the  Lloyd  Model  are  mathematical  related 
and  shown  to  be  equivalent  in  the  low  frequency  limit.  A  sample  first  generation  TIS  is 
used  to  compare  these  two  models  and  FLIR92  Model,  which  is  discussed  in  Chapter  HE. 
This  sample  first  generation  TIS  is  chosen  from  Reference  1,  presented  in  Appendix  A. 

B.  EYE-BRAIN  MODELS 

Several  researchers  have  tried  to  model  the  MTF  of  the  eye-brain  combination. 
The  Lloyd  and  the  Ratches  models  accept  the  eye-brain  combination  as  a  matched  filter. 

"  The  matched  filter  which  maximizes  the  signal-to-noise  ratio  (SNR)  (signal  being  the 
magnitude  of  the  output  from  the  matched  filter  and  noise  being  the  standard  deviation  of 
the  noise  fluctuations  )  at  a  time  q  for  the  case  that  the  noise  is  additive  (independent  of 
signal)  and  white  (the  power  spectrum  equals  a  constant  at  all  frequencies).  Note  that  for 
the  case  of  a  symmetrical  signal  and  for  ti  equal  to  zero,  the  matched  filter  has  precisely 
the  same  shape  as  the  signal.  (In  general,  the  matched  filter  is  the  mirror  image  of  the 
signal.)  "  [Ref.6,  p.l59]. 

FLIR92  uses  a  synchronous  integrator  model  for  MRTD  prediction  and  matched 
filter  concept  for  MDTD  prediction.  The  synchronous  integrator  model  can  be  defined  in 
the  following  way:  "  the  eye/brain  combination  will  integrate  over  the  entire  area  of  an 
image  even  though  the  image  has  been  smeared  out  over  a  large  distance  by  finite 
apertures.  In  the  formulation  for  signal,  the  integration  limits  are  plus  or  minus  infinity 
although  as  a  practical  matter  the  effective  distance  is  usually  much  smaller  because 
signal  integrated  in  the  low-amplitude  tails  of  a  blurred  image  in  a  real  imaging  system 
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increases  only  very  slowly  with  increase  in  integration  distance  from  the  image  center  or 
core  "  [Ref.  6,  p.  159].  The  matched  filter  concept  was  mentioned  in  Chapter  H. 

The  Visibility  Model  [Ref.  8]  uses  an  MTF^y^  based  on  an  approximate  sine  wave 
response  (SWR)  model.  This  eye  MTF  model  was  mentioned  in  Chapter  I  and  Figure  I.l 
shows  the  MTF  eye  model  as  a  band  pass  filter,  which  compares  well  with  cited  human 
eye  responses,  after  conversions  to  cycles/degree. 

C.  NETD 


Reference  6  derives  and  presents  the  NETD  in  the  following  way; 


NETD  = 


4fL(Afn)°-^ 

Tio(X)^D*(A.)dX 


[°C] 


(III) 


Afn  is  noise  equivalent  electrical  bandwidth,  is  the  area  of  a  single  detector,  Tio(?^p)  is 
the  optical  efficiency  of  the  lens,  and  Lx.  is  spectral  radiance  from  the  source,  f^^  is  the 
ratio  of  the  focal  length  of  the  collecting  lens  to  the  diameter  of  collecting  lens.  For 

calculation  by  hand,  is  given  5.2  10'®,  between  3.2  and  4.8)lm,  and  7.4  10 between 

8  and  13|lm.  Reference  1  derives  and  uses  some  approximations  to  give  a  simple  NETD 
equation,  which  is  shown  in  Equation  (II.2): 


NETD  = 


aPAo-CoDTXp)^ 


(11.2) 


In  this  equation  a  and  b  are  detector  dimensions,  area  of  the  optics.  To  is  optical 
transmission,a  and  p  are  detector  angular  subtenses,  is  1.48  lO  '*  between  8  and 
1 1.5|j.m. 
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D.  MRTD 

The  Ratches  MRTD  Model  and  the  Lloyd  MRTD  Model  define  the  MRTD  value 
in  the  following  ways: 


If  these  two  MRTD  models  are  compared  to  each  other,  it  is  seen  that  the  Ratches  model 

cc  7 

defined  and  —  =  L.  By  using  these  corresponding  terms  in  the  Equation  (11.3), 
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it  can  be  rewritten  as  follows: 


MRTDr 


^lSNR„NETDfT  /f  /f 


Hsys(fx)J^  H^ys(fy)H^(fy)dfy  V'^d  te  FAfr 


Si(lT) 


Substituting  related  terms  in  Equation  (11.6)  with  Equation  (II.4)  gives, 


MRTDr  =  MRIDtifT^-,  ' 

^  J-oc 


(fy)H^(fy)dfy 


[J-OO  Jo  '  J^eleCt(fx)H(Jj5p(fx,  fy)H^(fx)H'p(fy)Hsys(fy)dfx  dfy]®'^ 

Si  (It) 


Length  of  the  four  bar  pattern  is  seven  times  its  width, 


2  _  /2W 

7  V  L 


MRTDr  =  MRTDl 


JZ,  H?ys(fy)H2(fy)dfy' 


Sj  (it) 


Equation  (II.9)  shows  that  there  is  a  conversion  factor  between  the  MRTD„  and  MRTD 

R  L 
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which  shows  that  at  low  frequencies  these  two  models  are  equal  to  each  other.  Equation 
(11.4)  shows  as  spatial  frequency  goes  to  zero  MRTD;^  goes  to  zero,  so  MRTD,^  goes  to 
zero. 
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FLIR92  model  is  covered  in  the  next  chapter  and  Equation  (III.  19)  shows  that 
FLIR92  MRTD  goes  to  zero  when  spatial  frequency  goes  to  zero. 

The  Visibility  Model  agrees  with  the  Static  Performance  Model  except  at  very  low 
and  very  high  frequencies.  This  model  assumes  that  the  bar  pattern  acts  like  a  step 
function  in  the  zero  frequency  limit.  The  low  frequency  limit,  accepted  by  this  model, 
corresponds  to  a  minimum  system  degraded  temperature  difference,  which  is  enough  to 

be  recognized  by  observer  or  sensor.  This  model  proposes  that  there  is  a  critical  response, 
AT sc ,  below  which  the  four  bar  target  recognition  is  not  possible.  This  value  is  the 
response  amplitude  related  with  the  MRTD  of  the  measurement.  A  four  bar  frequency 

response  can  be  written  in  the  following  fcrm: 

“  MRTD  (n.16) 

By  using  Equation  (11. 16),  the  MRTD  zero  frequency  limit  relation  can  be  shown  as 
follows: 

MRTD(f  =  0)  =  jf^  (n.l7) 

In  this  part,  the  sample  system  which  is  given  in  Appendix  A  was  used  to  compare 
FLIR92  TIS  Performance  Model  to  the  Ratches  Model  and  the  Lloyd  Model  for  first 
generation  TISs.  Figure  II.  1.  shows  output  from  the  Ratches  Model  and  Lloyd  Model;  it 
can  be  seen  that  MRTD^  value  is  more  optimistic  than  MRTD^.  Figure  II.2.  compares  the 
Ratches  Model,  Lloyd  Model,  and  FLIR92  Model.  Figure  II.3  is  the  logarithmic 
presentation  of  Figure  II.2.  The  Ratches  Static  Performance  Model  was  developed  into 
I LIR90  and  finally  FLIR92.  These  figures  show  that  between  MRTD^  and  FLIR92 
MRTD  values  there  is  difference;  FLIR92  is  more  optimistic,  and  FLIR92  values  are 
closer  to  the  MRTDl. 

As  shown  the  Ratches  and  the  Lloyd  MRTD  Models  go  to  zero  when  spatial 
frequency  goes  to  zero  Since  FLIR92  MRTD  output  list  starts  from  0.05  cyc/mrad,  the 
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low  frequency  MRTD  limit  for  FLIR92  MRTD  can  not  be  established  from  the  data 
generated.  However,  Equation  (III.  19)  shows  it  goes  to  zero  for  zero  spatial  frequency, 
and  for  this  sample  system  it  determines  the  Nyquist  frequency  as  0.5  cyc/mrad,  which  is 
MRTD  values  are  not  calculated  after  this  frequency.  The  FLIR92  output  list  gives 
Horizontal  MRTD,  Vertical  MRTD,  and  2D  MRTD.  In  Figure  II.2  and  II.3  the  2D  MRTD 
was  plotted. 

MRTDr  and  MRTD^  are  shown  plotted  and  the  logical  flowchart  of  these  models 
is  given  in  Figure  II.4 
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Figure  II.  1  Ratches  and  Lloyd  MRTD  Models 
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Figure  II.4.  Visibility  Model  Logic  Flowchart 
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III.  U.S.  ARMY  NVSED  FLIR92  TIS 
PERFORMANCE  MODEL 


U.S.  Army  Night  Vision  Laboratory  developed  the  1975  Night  Vision  Laboratory 
Static  Performance  Model  for  first  generation  serial  and  parallel  scanning  TISs.  This 
model  was  one  dimensional  and  did  not  take  account  of  noise  sources,  except  for  random 
noise,  and  effects  of  sampling.  Because  of  the  improvement  in  TISs,  second  generation 
systems  replaced  first  generation  systems,  and  this  model  couldn't  predict  exactly  the 
second  generation  TISs'  performance.  In  1990,  NVEO  developed  the  FLIR90  TISs 
model  to  cover  the  different  noise  sources  and  sampling  effects.  This  model  was  updated 
as  the  FLIR92  TIS  Performance  Model  in  1992.  This  new  model  predicts  NETD,  MRTD, 
and  MDTD  [Ref.2]  [Ref.  9] 

This  model  can  be  defined  in  the  following  way;  "  FLIR92  models  parallel  scan, 
and  staring  thermal  imagers  that  operate  in  the  mid  and  far  infrared  spectral  bands.  The 
model  can  only  be  used  for  thermal  imagers  and  is  not  capable  of  predicting  performance 
for  any  other  class  of  electro-optical  sensor.  The  model  doesn't  predict  target  acquisition/ 
discrimination  range  performance."  [Ref.  9,  p.  ARG-1] 

In  this  chapter,  FLIR92  will  be  covered  under  three  headings: 

a.  MTFs 

b.  Noise 

c.  MRTD  and  MDTD 

A.  MTFs 

PLIR92  calculates  MTFs  in  three  main  groups.  The  product  of  these  groups  gives 
us  the  system's  overall  MTF.  MTFs,  which  are  calculated  by  this  model,  are  as  follows: 

1.  Prefilter  MTFs 

a.  Optics  MTFs 
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(1)  Diffraction-limited  MTF 

(2)  Geometric  Blur  MTF 

(3)  Measured  Optics  MTF 

b.  Detector  Spatial  MTF 

c.  Focal  Plane  Array  Integration  Time 

d.  Sample-scene  Phase  MTF 

e.  Image  Motion  MTFs 

(1)  Linear  Image  Motion  MTF 

(2)  Random  Image  Motion  MTF 

(3)  Sinusoidal  Image  Motion  MTF 

f.  Spare  Filter  MTF 

2.  Temporal  Postfilter  MTFs 

a.  Detector  Temporal  MTF 

b.  Electronics  Low  Frequency  Response 

c.  Electronics  High  Frequency  Response 

d.  Boosting  MTF 

e.  Spare  Temporal  Postfilter  MTF 

3.  Spatial  Postfilter  MTFs 

a.  Electro-optical  Multiplexor  MTF 

b.  Digital  Filter  MTF 

c.  Display  MTF 

(1)  CRT  Display  MTF 

d.  CCD  Charge  Transfer  Efficiency  MTF 

e.  Display  Sample  and  Hold  MTF 

f.  Eye  MTF 

(1)  Non-limiting  Eye  MTF 

(2)  Limiting  Eye  MTF 

g.  Spare  Spatial  Postfilter  MTFs 
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Equations  for  these  MTFs  are  given  in  Appendix  B.  The  short  listing  output  of  the 
FLIR92  gives  MTF  values  of  three  main  groups;  the  long  listing  output  of  FLIR92  gives 
all  MTF  values  under  the  three  main  groups. 

B.  NOISE 

NETD  was  used  to  predict  the  system  performance  for  first  generation  TISs. 

When  second  generation  TISs  began  to  be  used,  it  was  noticed  that  system  performance 
couldn't  exactly  be  predicted  by  using  NETD,  because  these  advanced  systems  use 
complex  signal  processing  techniques  which  cause  new  types  of  noises  and  these  noises 
can  be  defined  in  coordinate  system  by  averaging  them  in  specific  directions  [Ref.  14]. 
These  noises  affect  the  calculation  of  MRTD  using  the  performance  prediction  models 
which  were  used  for  first  generation  systems.  That  does  not  mean  that  first  generation 
systems  do  not  have  directional  noise,  but  detector  noise  is  the  dominant  noise,  affecting 
the  calculation  of  MRTD  for  these  systems. 

FLIR92  model  uses  the  3-D  noise  concept  for  scanning  and  staring  sensors. 
Dimensions  of  this  3-D  noise  are  t,  frame  number,  v,  row  number,  and  h,  column  number 
as  seen  in  Figure  HI.  1 . 


V,  row 

Figure  in.l.  3-D  Noise  Directions  (  After  Ref  [1 1] ) 
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There  are  totally  seven  noise  components  in  these  three  directions,  but  the  grand 
mean  of  all  these  seven  noise  components  can  be  treated  as  an  eighth  noise  component.  In 
the  following  Table  III.l,  these  eight  noise  components  are  shown. 


Table  III.l.  3-D  Noise  Components 
The  total  system  noise  U(t,v,h)  can  be  written  as  [Ref  9] 

U(f  v,h)=N^h  +N,^+Ntf,+Ny^+N^+Nh+N,+  s  (III.  1 ) 

where  all  terms  except  the  grand  mean  represent  fluctuations  about  a  mean  of  zero. 

1.  Directional  Averaging  Operators  (  D  Operators) 

3-D  Noise  components  are  random  fluctuating  values  about  the  mean  value  in  the 
direction  to  which  they  relate.  That's  why  we  need  to  average  these  values.  Variations  can 
be  integrated  out  by  using  D  operators  in  the  needed  directions.  Basically,  these  D 
operators  cancel  the  unwanted  noise  effect  components.  The  processing  definitions  of 
these  operators  can  be  seen  in  Table  IIL2. 


Noise  Term 

3-D  Process  Definition 

Ntvh 

[(l-DO(l-D,)(l-D^]{U(t,v,h)} 

Nvh 

[D,(l-D,)(l-DJ]{U(t,v,h)} 

Ntv 

[(l-DJ(l-D,)DJ{U(t,v,h)} 

N, 

[D,(l-D,)DJ{U(t,v,h)} 

Na. 

[(l-DOD,(l-D,J]{U(t,v,h)} 

Nh 

[DA(l-Dh)]{U(t,v,h)} 

N. 

[(l-DJDJDJ{U(t,v,h)} 

s 

[D.DJ)J{U(t,v,h)} 

Table  III.2.  Process  Definition  of  3-D  Operators  ( After  Ref.[12] ) 
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At  this  point  an  illustration  of  the  use  of  the  D  operators  is  presented.  If  N,^  is  to  be 
calculated,  (l-DJ,  (1-DJ,  and  D,,  are  used  since  we  want  to  cancel  the  noise  components 
only  in  the  h  (horizontal)  direction.  If  these  operators,  (l-D^ ),  (1-D^ ),  and  Dj,  are 
multiplied  with  the  system  total  noise  U(t,v,h),  it  will  give  the  desired  component. 
The  following  steps  are  provided  for  justification  of  the  process.  Which  is  illustrated  in 
Figure  ELS  [Ref  1 1] 

a.  The  D;  operator  has  0  value  for  operation  on  all  noises  in  the  temporal 
variation.  To  produce  the  (l-DJ  operator  it  is  necessary  to  subtract  all  noise  values  of  Dj 
from  1. 

b.  The  D^,  operator  has  0  value  for  operation  on  all  noises  in  the  vertical 
variation.  To  produce  (1-D^)  operator  it  is  necessary  to  subtract  all  noise  values  of 
from  1. 

c.  The  operator  has  0  value  for  operation  on  ail  the  noises  in  the 
horizontal  variation. 

d.  Multiply  the  same  noise  components  of  the  (1-D,),  (l-D^),  and  D^. 

e.  The  result  of  the  step-d  gives  only  noise  component  and  cancels  the 
other  noise  components.  Which  demonstrates  that  the  validity  of  the  row  of  Table 

ni.2. 

Another  approach  to  the  usage  of  D  operators  has  been  suggested  [Ref  12].  This 
approach  is  described  briefly  in  Appendix  B. 
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Figure  III.2.  The  Processing  of  D  Operators:  Calculation  of 
N,y,  the  Temporal  Vertical  Noise  Component 
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2.  Analysis 


Another  step  that  must  be  taken  in  the  3-D  noise  application  is  the  calculation  of 

the  standard  deviation  of  the  noise  components  expressed  in  degrees.  The  values  of  these 
standard  deviations  are  represented  by  <5^^^  ( where  xyz  is  the  related  direction ).  For 

example  calculation  of  the  Otvhis  shown  by  Equation  III.2. 

atvh  =  (stddev(N,vh  ( 1  -  F)))R  (in.2) 

where  the  average  standard  deviation  is  performed  in  this  case  over  directions  t,  v,  h 
along  pixel  streams.  In  this  equation,  F  is  the  low  frequency  trend,  which  is  low  frequency 
fixed  pattern  noise  [Ref.  14],  and  R  is  the  ratio  of  input  temperature  to  the  output 
intensity  measured  in  (“CA^olt).  High  frequency  variations  have  the  biggest  effect  on  the 
system  performance  predictions  so  it  is  necessary  to  remove  the  low  frequency  trend  in 
spatial  noise.  By  removing  the  low  frequency  trend,  high  frequency  variations  will  be 
made  dominant  in  the  direction  of  the  noise  component.  The  values  of  the  other  noise 
components  can  be  calculated  by  changing  Equation  in.2  for  their  directions.  For 
calculation  of  the  s-component,  we  don’t  need  any  trend  subtraction  or  standard  deviation; 
we  multiply  it  by  ratio  of  input  temperature  and  output  intensity.  Table  HI.S  shows  3-D 
noise  components  and  their  sources  [Ref.  2,  p.386].  The  following  equation  shows  the 
relation  between  the  classic  NETD  and  3-D  noise  component  [Ref.  9]; 


atvh  =  NETD  X 


4  fro 


(ni.3) 


which  follows  from  the  fact  that  Gtvh  is  a  composite  measure  of  noise  for  all  three 

directions.  Afp  is  actual  system  noise  bandwidth  associated  with  the  system  electronics 
prior  to  display;  for  scanning  and  staring  systems  this  bandwidth  can  be  shown  in  the 

following  way: 
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3D  Noise 
Component 

Description 

Serial  Scan 

Parallel  Scan 

Staring  Array 

CT/v/a 

Random  3-D 

noise 

Random  and  1/f 
noise 

Random  and  1/f 
noise 

Random 

Spatial  noise  that 
does  not  change 
from 

frame-to-frame 

- 

- 

FPN 

^th 

Variations  in 

column  averages 
that  change  from 
frame-to-frame 
(rain) 

Microphonics 

Microphonics 

Readout  noise 

(5  tv 

Variations  in 

row  averages 

that  change  from 
frame-to-frame 
(streaking) 

1/f  noise 

Transients 
(flashing 
detectors),  1/f 
noise 

Readout  noise 

CTv 

Variations  in 

row  averages 

that  are  fixed  in 
time  (horizontal 
lines  or  bands) 

Line-to-line 

interpolation 

Detector 

gain/level 

variations, 

line-to-line 

interpolation 

Readout  noise, 

line-to-line 

interpolation 

Oh 

Variations  in 

column  averages 
that  are  fixed  in 
time  (vertical 

lines) 

Shading 

Shading 

Readout  noise 

Ot 

Frame-to-frame 

intensity 

variations 

(flicker) 

Frame 

processing 

Frame 

processing 

Frame 

processing 

Table  III.3.  Seven  Noise  Components  of  The  3-D  Noise  Model  [Reference  2] 
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for  scanning  systems. 


Afp  =  J“  S(v)H?pF(v)dv  (in.4) 

for  staring  systems, 

Afp  =  j:s(v)(21?i)dv  (1115) 

where  t;  is  the  focal  plane  array  integration  time,  S(v)  is  normalized  detector  noise  power 
spectrum  and  H^-pp  is  the  temporal  post  filter  defined  in  Appendix  B.  Af„  is  the 
equivalent  noise  bandwidth  for  the  NETD  and  can  be  defined  as: 

Af„  =  J“  S(v)H?,f(v)dv  (in.6) 

where  is  the  standard  NETD  reference  filter  [Ref.  1]. 

The  total  system  noise,  generated  from  the  statistically  independent  standard 
deviation  noise  components,  <5^^^ ,  can  be  written  as  [Ref.  13,  p.  ARG-33]: 

^  +  Cfvh  +  Cfh  +  (in.7) 

Total  system  noise  in  horizontal  and  vertical  directions  can  be  seen  in  the  following 
Equations  (III.8)  and  (111.9): 

for  the  horizontal  direction, 

Qh  =  [af,hEtEv(f)Eh(f)  +  a,\Ev(f)Eh(f)  +  a^E,Eh(f)  +  a^Eh(f)]°'  (ffl.S) 
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for  the  vertical  direction. 


=  [afvhEtEv(f)Eh(f)  +  avhEv(f)Eh(f)  +  G^vEtEv(f)  +  a?Ev(f)]°'^  (111.9) 

Ej,  E^(f),  and  E[,(f)  are  the  integration  factors  of  the  display/eye/brain  in  the  related 
directions.  These  integration  factors  are  [Ref.  13]: 


The  temporal  integration  factor  is. 


Et- 


CCt 

Fr  te 


(in.  10) 


where  is  the  system  frame  rate,  te  is  the  eye  integration  time,  and  ttt  is  the  sample 
correlation  factor.  CC^  is  usually  taken  to  be  1 . 


The  vertical  integration  factor  is. 


Ev(f) 


gy 

Rv  Lv(f) 


(in.ii) 


where  is  the  vertical  sampling  rate  (samp./mrad),  L^(f)  is  the  vertical  spatial  integration 
limit,  which  in  this  approximate  form  is  the  vertical  dimension  of  the  MRTD  bar  target, 
and  ttv  is  the  vertical  sample  correlation  factor  [Ref.  9]. 

The  horizontal  integration  factor  is. 


Eh(f)- 


«h 

RhLh(f) 


(in.  12) 


where  R^,  is  the  horizontal  sampling  rate  (samp./mrad),  Lh(f)  is  the  horizontal  spatial 

integration  limit,  which  in  this  approximation  is  the  horizontal  dimension  of  the  MRTD 
bar  target,  and  tth  is  horizontal  sample  correlation  factor.  For  staring  systems  (Xv  and  Gh 
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are  usually  taken  to  be  one  [Ref.  9]. 


3.  Noise-Data  Input  Groups 

There  are  three  different  noise-data  input  groups  in  FLIR92.  By  using  any  one  of 
these  groups,  desired  noise  information  can  be  entered  into  the  model.  These  noise-data 
groups  are: 

a.  3-D  Noise  Default:  If  we  don't  specifically  assign  noise  values,  FLIR92 
model  assigns  default  values  as  a  percentage  of  atvh  for  the  system's  critical  noise 
components,  which  are  shown  in  Table  III.4.  Moderate  noise  level,  for  scanning  systems, 

and  low  noise  level,  for  staring  systems,  are  recommended  and  used  as  default  values. 

b.  3-D  Noise  Measurements:  This  data  group  is  used  to  enter  the  exact 
3-D  noise  values. 

c.  3-D  Noise  Estimates:  This  data  group  is  used  to  enter  the  estimated  3-D 
noise  values  as  a  percentage  of  atvh  • 


System 

Noise 

Component 

No  Noise 
Default 

Low  Noise 
Default 

Moderate 

Noise 

Default 

High  Noise 
Default 

Scanning 

0 

0.25  X  atvh 

0.75  X  atvh 

1 .0  X  O' tvh 

Staring 

CJvh 

0 

0.40  X  atvh 

- 

- 

Table  in.4.  FLIR92  3-D  Noise  Default  Values 


C.  MRTD  AND  MDTD 


FLIR92  uses  synchronous  integrator  model  for  MRTD  prediction  and  matched 
filter  concept  for  MDTD  prediction.  Synchronous  integrator  model  and  matched  filter 
concept  were  mentioned  in  Chapter  II. 

FLIR92  corresponding  to  the  direction  of  the  standard  four  bar  pattern  calculates 
the  horizontal  and  the  vertical  MRTD.  Horizontal  and  vertical  orientations  of  the  four  bar 
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pattern  are  shown  in  Figure  111.3.  Using  these  two  MRTD  values,  FLIR92  calculates  the 
MRTD^d,  for  which  the  spatial  frequency  is  defined  by  Equation  (IE.  13). 


flD  =  Jfx  ■  fy 


(in.  13) 


- ►  h  ' - ►  h 

a.  horizontal  b.  vertical 

Figure  in.3.  MRTD  Target  Patterns  (From  Ref.  9,  p.  ARG-2) 


MRTD  and  MDTD  are  described  by  the  following  equations  [Ref.  9]: 


MRTDz(f)  = 


Hsysz(f) 


(ni.i4) 


where  the  z  subscript  refers  to  the  horizontal  or  vertical  direction  as  required  and  MDTD 

is 


MDTD(f)  = 


SNRthOtvh^ 

AiQhtfiQvCf) 


(in.  15) 


We  can  write  the  total  system  noise,  Equation  (111.8)  and  (III.9),  in  the  form 


^^z  =  atvhkz  [EtEh(f)Ev(f)]°-' 


(ni.i6) 
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where  is  the  correction  factor  in  the  horizontal  or  vertical  direction.  Using  Equations 
(III.  14),  (III.  15),  and  (III.  16),  MRTD  and  MDTD  can  be  rewritten  as  follows: 


MRTD,(f)  =  [£f  (ni.i7) 

If  horizontal  MRTD  for  scanning  systems  is  written  by  using  the  exact  equations  of  the 
display/eye^rain  integration,  [Ref  9,p.ARG-39]  instead  of  using  any  approximation,  it 
gives  the  following  equation; 


2  f 

J  5^(v)«„y{v) 


•  /'tiv 

2 

sm(^ 

1 

dv 

TtV 

\  2/ 

) 

J 

(in.i8) 


X 


7  71  V 

\  Tr  J 


2  ^ 

dv 


10.5 


By  choosing  one  of  the  integrals  in  (HI.  18),  i.e.  first  integral,  it  can  be  shown  that 
FLIR92  MRTD  value  goes  to  zero.  If  ^  is  defined  as  x,  the  first  integral  can  be  written 
as  follows: 


^  ]  S(v)H„Ay) 


(in.i9) 


Equation  (HI.  19)  proves  that  FLIR92  MRTD  is  equal  to  zero  when  spatial  frequency  is 
equal  to  zero.  MDTD  value  for  FLIR92  is, 

MDTmf)  =  [E,E^(f)Erm'‘''  (ni-20) 

where  Qh(f)  and  Qv(f)  are  given  in  following  equation; 
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Q»(f)  =  J  H|,,(v) 


2 

dv 


(in.2i) 


By  using  Equation  (111.8),  (in.9),  and  (III.  16)  correction  factor,  k^(f),  can  be  written  : 
for  the  horizontal  direction  as 


(in.22) 


for  vertical  direction  as 


kv(f) 


(in.23) 


If  it  is  assumed  that  the  recommended  default  values,  i.e.  moderate  noise  for 
scanning  systems  and  low  noise  for  staring  systems,  are  used  from  Table  111.4,  then  the 
correction  factors  can  be  simplified  to  the  following  forms: 
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for  scanning  systems. 


kh(f)  =  1  (in.25) 


k.(f)  =  kMDTD(f)  =  [l  +  (0.5625)^  +  (0.5625)^f (m.26) 

for  staring  systems, 

kh(f)  =  kv(f)  =  kMDTD(f)  =  [  1  +  (0.16)^]°'^  (m.27) 

These  correction  factors  can  be  inserted  in  to  the  MRTD  and  MDTD  equations,  (HI.  17) 
andCffl.lS). 
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IV.  ALIASING 


Analog-to-digital  and  digital-to-analog  signal  conversions  cause  degradation  in 
TISs.  Aliasing  and  blurring  contribute  significantly  to  this  degradation.  Aliasing  is 
created  by  sampling  at  less  than  the  Nyquist  frequency.  Blurring,  a  term  commonly  used 
for  images,  is  generated  when  the  image  is  spatially  low  pass  filtered  [Ref  38].  In  this 
chapter,  aliasing  models  relevant  to  electronic  imaging  will  be  covered.  Figure  IV.l 
shows  a  block  diagram  of  a  sampled  TIS.  This  model  will  be  described  in  more  detail  in 
the  next  section. 


Figure  IV.l.  Block  Diagram  of  a  Sampled  TIS  ( After  Ref  [17] ) 


A.  MODEL 


The  processes  illustrated  in  Figure  IV.  1  are  described  in  the  space  domain  and 
frequency  domain  by  following  steps  [  Ref.  17  ] : 


41 


d.  The  reconstruction  process  is  the  convolution  of  the  sampled  image 
signal,  Ij ,  with  the  transfer  function  of  the  reconstruction  process,  h^  (x,y): 

I  (x,y)  =  Is  (x,y;  X,Y)*h,  (x,y)  (IV.3) 

2.  In  the  frequency  domain,  these  steps  (a,  b,  c,  d)  can  be  expressed  by  using 
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Fourier  transforms.  The  two  dimensional  Fourier  transform  operator  can  be  defined 
[Ref.38]: 

oo  oo 

ci2'{G(x,y)}  =  J  J  =  G(fx,fy)  (rv.4) 

—oo  — oo 

where  4  and  f^  are  the  spatial  frequencies.  The  inverse  transform  is  at  the  same 
mathematical  form  except  the  complex  exponentials  are  complex  conjugated.  This 
relation  can  also  define  the  one  dimensional  transform.  If  G(x,y)  separates,  i.e. 

G(x,y)=  G(x)G(y)  then  c52{G(x,y)}=  ^d{G(x)}  62',^  {G(y)}. 

a.  The  Fourier  transform  of  the  scene  is: 

r(fx,fy)=e$?{I(x,y)}  (IV.5) 

b.  The  image  spectmm  is  obtained  by  taking  the  Fourier  transform  of  the 
image  signal  I^,(x,y).  Application  of  the  convolution  theorem  [Ref.  38]  leads  to: 

U(fx,fy)  =  I  (fx,fy)Hi(fx,fy)  (IV.6) 

c.  The  spectrum  of  the  ideal  sampling  waveform  [Ref.  17]  is: 

Ps(Xfx,Yfy)  =  3^^^^J^8  (fx-f  j,(^fy-f  j  (IV.7) 

where  m  and  n  are  indices  that  identify  where  sampling  occurs.lt  is  worth  noting  that  the 

ideal  case  can  be  applied  to  nonideal  sampling,  i.e.  sampling  with  pulses  instead  of 
impulses,  by  prefiltering  the  image  spectrum  I  (fx ,  fy )  with  a  lowpass  transfer  function 
associated  with  the  window  process  of  the  detector  [Ref.  1.  p.  376]. 
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d.  The  ideal  sampled  image  is. 


I  s(fx, fy ; X,  Y)  =  I  fy)  *  Ps(Xfx,  Yfy)  (IV.8) 

By  using  Equation  (IV.7),  the  sampled  image  spectrum  can  be  expressed  as  an  infinite 
sum  of  terms  as  shown. 

I  S(f.,  f,;  X.  Y)  =  .[(fx  -  f  ),  (f,  -  f )]  (IV.9a) 

or  equivalently 

= T  m(fx ,  f ,) + J__  (fx  -  f  ).  (^f ,  -  f  ) j  (iv,9b) 

(m,  n)  0,  0) 

where  the  factor  (XY)''  common  to  all  terms  has  been  dropped.  Equation  (IV.9a)  and 
(IV.9b)  can  be  presented  as  follows: 

I  s(fx,  fy ;  X,  Y)  =T  m(fx,  fy)  +T  ma(fx,  fy ;  X,  Y)  (IV.  10) 

where  the  addition  of  the  subscript  "a"  is  intended  to  indicate  aliasing.  After  substituting 
(IV .  10)  back  into  (IV.  3)  and  taking  the  Fourier  transform  leads  to  the  reconstructed 
image  spectrum: 

Ir(fx,fy)=Tn,(fx,fy)Hr(fx,fy)  +  Ima(fx,fy;X,Y)Hr(fx,fy)  (IV.ll) 
The  inverse  fourier  transform  of  the  Equation  (IV.  11)  gives  the  reconstructed  image 


signal  in  the  space  domain,  which  is  equal  to  Equation  (IV.S): 


;  (x,y)  =  (x,y)  + 1^3,  (x,y)  (IV.12) 

(x,y)  is  the  ideal  reconstructed  image  without  aliasing  and  the 

reconstruction  of  the  higher  order  terms  in  which  occurs  due  to  sampling  [Ref.  38].  The 
aliasing  effect  due  to  the  addition  of  see  Equation  (IV.  12),  can  be  considered  as  noise 
[Ref.  17].  Substitution  of  (IV.9b)  into  (IV.  10)  followed  by  (IV.  11)  shows  that 

I  „(fx,  fy)  =T  „(fx.  fy)  H,(fx.f,)  av.l3) 


or  from  (IV.6) 


I  „.(fx,fy)=  I(fx,fy)Hi(fx,fy)Hr(fx,fy) 


(IV.  14) 


which  can  be  referred  as  the  DC  term  in  (IV.12)  since  no  frequency  shifting  is  involved. 
The  substitution  process  described  above  also  demonstrates  that 


I™(fx,fy)=  I  X  I„(fx-f,fy-7)Hi(fx-f,fy-f)H,{fx.fy) 


m=-oo  n=-oo 


(IV.15) 


which  includes  an  infinite  number  of  terms.  From  this  point  only  the  terms  with 
|n|  =  |m|  =  1  will  be  considered  significant.  Higher  order  terms,  i.e.  n  >  1  m  >1, 
appear  at  successively  higher  multiples  of  the  sampling  frequencies  (  1/X,  1/Y).  Since 

these  terms  are  multiplied  in  (IV.  15)  by  H^f^,  fy),  which  is  known  to  be  a  low  pass  filter, 
they  are  neglected  in  the  subsequent  analysis. 


45 


B.  ANALYSIS  AND  SIMULATION  OF  THE  ALIASING 


The  thermal  image  system,  denoted  as  sample  system  A,  has  parameters  which  are 
defined  in  Appendix  A.  The  specifications  for  this  system  have  been  used  to  generate 

both  the  DC  reconstruction  transfer  function  H/f^,fj,),  see  Equation  (IV.  14),  and 

the  first  order  |n|  =  |m|  =  1  transfer  function  from  (IV.  15).  As  previously  mentioned  the 

higher  orders,  i.e  |n|  >  1  |m|  >  1 ,  are  neglected.  To  demonstrate  the  aliasing  effect  the 
equations  described  in  the  previous  section  will  be  applied  to  the  parallel  scanning  TIS, 

system  A,  for  which  only  one  spatial  frequency,  4,  is  significant.  An  adjustment  in 
notation  is  introduced  according  to: 

Hpre(fx)  =  Hi  =  Hopt(fx)Hdet(fx)  (IV.16) 

and 


Hpost(fx)  —  Hr  —  Helect(fx)  H(jisp(fx)  (IV.  17) 

which  is  motivated  by  noting  that  previous  to  sampling  the  image  is  affected  by  both  the 
optics  and  the  averaging  effect  of  the  detector  window  [Ref.  1].  Also  after  sampling  the 
image  signal  is  affected  by  both  the  electronics  and  the  display.  It  follows  from  (IV.  15) 
that  the  first  order  term,  acting  as  noise,  appears  as  an  overlap  involving  (now  known 
as  Hp^j,)  and  a  shifted  version  of  Hj  (now  known  as  Hp^J.  This  motivates  the  definition: 

Hover  ~  Hpost(fx)  Hpre(fx  i  fs)  (IV.IS) 

where  the  ±  correspond  to  n=  ±  1 ,  in  (IV.  15),  respectively  and  f.=  1/X.  This  overlap 
transfer  function  will  be  applied  to  the  image  spectrum  I  ni(fx)  shifted  by  ±fs .  Figure  IV.2 
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shows  sketch  for  the  expected  frequency  dependence  of  magnitude  transfer  function  for 
Hpre,  see  Equation  (IV.  16),  see  Equation  (IV.  17),  and  Hover,  see  Equation  (IV.  18). 
Figure  IV.3  the  corresponding  curves  for  the  TIS  sample  system  A. 

1.  Sample  Scene  Phasing 

The  effect  of  spatially  shifting  the  sampling  function  relative  to  a  fixed  image  has 
been  referred  to  as  sample  scene  phasing.  For  purposes  of  evaluating  the  this  effect  on  the 
reconstructed  image  it  is  necessary  to  represent  the  ideal  sampling  waveform  and  its 
spectrum.  Figure  IV.4  defines  the  two  different  cases  of  interest,  i.e  ideal  sampling 
waveforms  with  and  without  spatial  shift.In  Case- A,  the  reference,  there  is  no  phase  shift 
due  to  in  the  sampling  process,  but  Case-B  has  a  phase  shift  due  to  a  spatial  shift  in  the 
sampling  process.  The  exponential  Fourier  series  coefficients,  F„,  can  be  obtained  for 
Case-B  by  first  analyzing  Case-A  and  then  using  the  shift  theorem  to  relate  spectrum  B  to 
spectrum  A.  For  any  periodic  function  [Ref  37]. 

¥  =  ^^{g(x)}/X  (IV.19) 

where  g(x)  is  1  cycle  of  the  periodic  waveform  centered  on  x=0  and  X  is  the  period  of  the 
waveform.  For  Case-A,  see  Figure  IV.4,  g(x)=  5(x)  and  X=l/fs.  It  follows  that 

Fi  =  f,  (rv.20) 

and  therefore  after  reference  to  Figure  IV.4  and  using  the  shift  theorem  [Ref  37] 

Fl  =  Fi  (IV.21) 

where  dx  is  the  spatial  shift  in  the  sampling  process.  Note  that  dx  can  be  limited  to 
0<dx<\lfs 
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Spatial  Frequency 


Figure  rv.2.  Presentation  of  Aliasing  by  Using  MTF 

a.  Presampling  and  Postsampling  Transfer  Function 

b.  Presampling  Function  with  Shifted  Replica 

c.  Alias  MTF  and  Ideal  MTF 


Figure  IV.3.  Hp^st(4),  Hp^^(f^-f^),  for  the  Sample  System  A 
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CASE-A 


CASE-B 


Figure  IV.4.  Two  Different  Cases  for  Sampling  Points 


Equation  (IV.  18)  can  be  separated  into  its  real  and  imaginary  parts: 

The  real  part  is  an  even  function, 

Re{F®  }  =  f s  cos(27rnfs  dx)  (IV.22) 

The  imaginary  part  is  an  odd  function, 

Im{Fn}  =  fs  sin(27infs  dx)  (IV.23) 

The  exponential  Fourier  series  representation  for  the  periodic  sampling  waveform  B  is 

fp(x)=  X  £  Fn  [cos(n27ifsx) +jsin(n27cfsx)]  (IV.24) 

^  n=— oo  n=— oo 

By  defining  F„  in  two  parts  as  real  and  imaginary.  Equations  (IV.22),  (rv.23).  Equation 
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(IV.  24)  can  be  written  in  the  following  form; 


f?(x)  =  [(Re{Fn}cos(27tnfsx)  -  Im{Fn}sin(27infsX))+  (IV.25) 

j(Im{Fn  }cos(27i;nfsX)  +  Re{Fn  }sin(27rnf  sx))] 

By  substituting  the  real  and  imaginary  F„  from  Equation  (IV.22)  and  (IV.23)  to  (IV.25) 
leads  to  a  simpler  expression 


f®(x)  =  fs  +  2  I  [cos(27infsdx)cos(27rnfsX)-sin(27rnfsdx)sin(27rnfsx)]  (IV.26) 

n=l 

where  a  factor  of  2  accounts  for  terms  having  a  negative  index  in  (IV.25).  When  n  is 
equal  to  zero,  (IV.25)  gives  the  DC  component  of  the  sampling  function: 


f?oW  =  fs  (IV.27) 

When  n  is  equal  to  one,  Equation  (rv.26)  gives  the  term,  which  is  used  in  the  subsequent 
calculations  to  estimate  the  aliasing  effect: 

f (x)  =  2  fs  [cos(27ifsdx)  cos(27tfsX)  -  sin(27ifsdx)sin(27ifsx)]  (IV.28) 

which  can  be  simplified  in  the  following  way; 

f^,  (x)  =  2 fs  sin[27tfs(x  +  dx)  +  f ]  (IV.29) 
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2.  Calculation  Of  The  Aliasing 


The  input  image  signal,  a  periodic  square  wave  of  period  fx^  is  given  in  Equation 
(IV.27)  [Ref. 8]: 

I(f,x)  =  ^+  i  i  sin(27tnfxx)  (IV.30) 

which  provides  for  MRTD  analysis,  a  convenient  mathematical  approximation  to  a 
four-bar  pattern.  Multiplication  of  the  first  order  term  from  the  sampling  function  which 
is  Equation  (IV.29),  with  the  input  image  signal,  Equation  (IV.30),  gives  the  first  order 
spectral  replica  of  the  image  signal,  which  is  centered  due  to  sampling.  Equation  (IV.27) 
shows  that  when  n  is  equal  to  zero,  the  DC  component  is  equal  to  f^ .  Because  the  DC 
spectral  replica  is  the  ideal  reference  which  is  not  included  in  the  following  equation,  this 
first  factor  is  non-essential,  and  we  can  write 

I''(x,  f)  =  sin(27cfs(x  -I-  dx)  -I-  f )  -f  2  S  H[cos(27i(fs  -  nf)  x  +  2jifsdx  +  f )  - 

^  n=l  ^ 

cos(27i(fs  +  nf)x  -h  27tfsdx  +  f )]  (IV.3 1) 

Equation  (IV.31)  can  be  simplified  by  defining  the  three  terms; 

Terml=sin(27rfs(x  +  dx)-i-|)  (IV.32a) 

oo 

Term2=2  Z  ^[cos(27T:(fs  -  nf)  x  +  27cfsdx  +  f )  (IV.32b) 

n=l,3,5...  ^ 

Term3=2  Z  i  [cos(27t:(fs  +  nf)x  +  27lfsdx  +  f) 

n=l,3,5...  ^ 


(IV.32C) 


It  can  then  be  is  rewritten  by  using  these  three  terms,  as 


I(x,f)=  Terml  +  Term2  -  TermS  (IV.32d) 

Since  the  image  is  filtered  by  Hp^  prior  to  sampling  the  alias  OTF,  shown  in  Figure  IV.2, 
needs  to  operate  on  the  shifted  image.  Figure  IV.5  defines  the  effects  of  a  transfer 
function,  H,  on  an  input  sinusoid. 


OTF  =  |Hpost(fx)|  |Hpre(fx  “  fs)+ZPTFpost(fx)]  (IV.33) 


Cos(2  T  a!  X  +  b  ) 

_ 

1  H(ai 

1  Cos(2¥  a' X  +  b’ +  PTF(a) ) 

- 

Figure  IV.5.  The  H-Rule  (a'  and  b'  are  arbitrary  constants) 


Equation  (IV.32.a-d)  and  (IV.33)  gives  the  aliasing  term,  which  can  be  written  by  using 
Terml,  Term2,  and  Term3. 

A  =  OTF'(Terml )  +  OTF'(Term2)  +  OTF'(Term3)  (IV.34) 


when 


OTF'(Terml)  =  ^|H2(fs)||  Hi(0)|  sin[27tfs(x  +  dx)  +  f +  PTF2(fs)  +  PTF,(0)] 

(rv.35a) 

OTF'(Term2)  =  1^  ^ lH2(fs  - nf)| |Hi (-nf)|cos(27i(fs  - nf)x  +  27tfsdx  +  f  (IV.35b) 

+PTF2(fs  -  nf)  +  PTFi  (-nf); 
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and 


0TF'(Term3)  =-  E  H|H2(fs  +  nf)|  |Hi (nf)lcos(27c(fs  +  nf)x  +  2jifsdx  +  f  (IV.35c) 

+PTF2(fs  +  nf)  +  PTF,(nf)) 

where  for  compactness  subscript  2  denotes  pre  and  subscript  1  denotes  post.  It  is 
important  to  note  that  consistent  with  the  symmetry  conditions  of  transfer  function 
[Ref.37]  the  magnitude  is  an  even  function  about  DC  and  the  phase  is  an  odd  function 
about  DC.  Equation  (IV.34)  shows  the  first  order  aliasing  term.  Higher  order  aliasing 
terms  are  not  significant  as  it  can  be  seen  from  Figure  IV.5.  In  Equations  (IV.35.a-c),  Hj 
and  Hj  are  used  for  the  terms  which  are  shown  in  the  following  equations.  PTEj  is  used 
for  detector  and  PTF2  is  used  for  electronics. 

H,  =  H,.,  (IV.3«a) 

and 

U,=  H,,,  (IV.36b) 

3.  Simulation  of  the  Aliasing 

The  sample  system,  which  is  given  in  Appendix  A,  is  used  to  simulate  the 
aliasing.  In  this  simulation  dx  is  0.1  mrad,  and  spatial  frequency  is  0.15  cyc/mrad.  Figure 
IV.5  shows  the  square  wave  image  signal  and  expected  image  output  signal  with  phase 
shift,  because  of  the  phase  transfer  function  of  the  system.  Figure  IV.6  presents  the 
expected  image  output  signal  and  aliasing,  which  was  defined  by  Equation  (IV.24). 

Figure  rv.7  shows  the  square  wave  image  signal  and  output  image  with  aliasing.  Figure 
IV. 8  is  the  presentation  of  the  square  wave  input  image  signal,  expected  image  output, 
and  output  image  with  aliasing.  These  figures  show  that  aliasing  can  be  considered  as  an 
additional  noise  and  causes  distortion  on  the  expected  output  image  signal. 
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normalized  response 


Square  Wave  and  Square  Wave  Response 


Figure  IV. 5.  Square  Wave  Image  Signal  and  Expected  Image  Output  Signal 

(After Ref.  [8]) 


Expected  Output  and  Aliasing  Effect 


Figure  IV.6.  Expected  Image  Output  Signal  and  Aliasing 
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normalized  response  m  S  normalized  response 


Square  Wave  and  Square  Wave  Response  w/Aliasing 


re  IV.7.  Square  Wave  Image  Signal  and  Image  Output  with  Aliasin 


iquare  Wave  and  Square  Wave  Response  w/o  and  w/Aliasing 


Figure  IV.8.  Square  Wave  Image  Signal,  Expected  Image  Output 


55 


C.  REDUCING  THE  ALIASING  EFFECT 


There  are  different  methods  to  reduce  the  amount  of  the  aliasing  ( Ref.  [18] ); 

1.  Increasing  the  fill  factor  to  100% 

2.  Increasing  the  sampling  rate  by  decreasing  the  FOV 

3.  Increasing  the  number  of  detectors 

4.  Microscanning 

The  first  three  of  these  methods  have  some  fabrication  problems  and  they  are  not 
preferred.  That  is  why  microscanning  will  be  presented  as  a  method  to  reduce  the  aliasing 
for  staring  systems  (  Ref.  [18] ). 

1.  Microscanning  Process 

In  this  technique,  FOV  of  the  detector  is  shifted  around  by  increments  relative  to 
the  scene.  This  shifting  distance  is  related  to  the  spacing  between  the  detectors,  and 
amount  of  the  shifting  depends  on  the  desired  microscanning  level.  Figure  IV.9  shows  the 
shifting  of  the  FOV. 


(a)  (b)  (c) 

Figure  IV.9.  Microscanning  Process 

Equation  (IV.2)  will  be  used  to  represent  the  sampled  image,  and  it  is  rewritten  in  the 
following  equation: 
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1,  ( x,y;  X,  Y )  =  U  x,y  )Ps[  (f  ) » (|) 

The  ideal  sampling  function  can  be  defined  in  the  following  way: 


(IV.37) 


(IV.38) 


where  the  comb  function  is  used  to  represent  an  infinite  array  of  uniformly  spaced 
impulses  [Ref.  38],  By  using  Equation  (IV.37)  and  (IV.38),  the  sampled  image  signal  can 
be  rewritten  as  follows: 


I.  ( ’‘.y;  X,Y  )  >  [^  b{f  ,  f)  *  Kx,y)  .  h,(x,y)_  x 


(IV.39) 


where  Equation  IV.  1  has  been  used  to  relate  the  image  to  the  scene.  In  this  equation 
is  the  shape  of  the  each  pixel  and  d^ ,  dy  are  the  widths  of  the  pixels  in  the 

\Ux  dy  J 

relative  directions.  If  four  microscan  steps,  shown  in  Figure  (IV.8),  are  applied  to  this 
image  signal,  the  linear  superposition  of  "snapshots"  taken  during  1  cycle  of  the 
microscan  process  can  be  represented  in  the  following  way: 


4(x,3;;X,  Y) 


1 

4 


^{combQ,  |)  +  comb(^-  |)+ 
comb{^,  f  -  5)  +  comb(^  ~2’y~  2^^ 


(rv.40) 


Equation  (IV.40)  shows  that  integration  time  is  increased  by  a  factor  of  4,  because  the 
image  signal  is  integrated  four  times  for  one  sampled  image  signal.  The  Fourier  transform 
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equation  is  written  as  follows: 


1  Y)  =  '{b(f,dyjyd,)~l  (fy,fy)Ht(fy,fy)  X  ;  (IV.41) 

comblfJC,fyY){  1  +  } 

By  using  Equation  (IV.  5),  it  can  be  defined  in  the  following  form: 

I  s(f„fy-,X,Y)=  \[bJfydy,fydy)~l(f„fy)H,(fy,f,)]  X  ;  {IV.42) 

m-—<x)  «=-<50  V  i  ) 

which  follows  after  noting  that  the  localized  contribution  of  the  impulses  permits  the  sum 
of  complex  exponentials  in  (IV.41)  evaluated  at  f=  j  and  fy=  j .  The  final  step  is 
Equation  (IV.42),  which  reduces  to  zero,  if  both  m  and  n  are  not  even.  It  removes  the  odd 

harmonics  and  decreases  the  amount  of  the  overlap.  Figure  (IV.9)  shows  this  process. 


a.  Harmonics  overlap  each  other  b.  Odd  Harmonics  are  removed 

Figure  IV.  10.  Removing  Odd  Harmonics  by  Microscanning 
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V.  LABORATORY  MEASUREMENTS 


To  compare  the  FLIR92  TISs  Performance  Model  to  the  measured  data  for 
second  generation  TISs,  two  laboratory  measurements  were  done.  An  Amber 
Engineering-  AE4128  IR  Imaging  System  and  a  Mitsubishi  Electronics  IR-M500  Imager 
were  used  for  these  experiments.  These  two  systems'  parameters  are  given  in  Appendix 
A.  Comparison  of  these  measurements  to  FLIR92  is  shown  in  Chapter  VI. 

A.  LABORATORY  SETUP 


Figure  V.  1  shows  the  laboratory  setup  used  for  both  of  the  measurements.  The 
oscilloscope,  which  is  shown  in  this  figure,  was  used  only  for  the  Mitsubishi  IR-M500 
experiment  as  an  Automatic  Target  Recognition  (ATR)  unit. 


Figure  V.l.  Laboratory  Setup 
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Aluminum  was  used  for  the  target  plates,  painted  in  a  non-reflective  flat  black  to 
have  uniform  emissivity.  A  variable  heater  was  used  to  heat  the  back  plate.  Front  plate 
was  placed  10  cm.  in  front  of  the  back  plate  to  keep  it  at  ambient  temperature.  Four  bar 
patterns  which  were  cut  through  the  front  plate  had  standard  size  with ,  length  equal  to 
seven  times  the  width,  to  give  different  spatial  frequencies  related  to  the  distance  from 
the  target  to  the  imager.  The  thermocouple  read-out  displays  had  one  digit  decimal 

precision;  to  have  more  precise  reading,  these  displays  were  used  in  Fahrenheit  mode, 
provided  ±0.0556°C  tolerance.  This  tolerance  limited  the  precision  of  the  measurements. 
The  target  apparatus  was  placed  in  front  of  a  non-reflective  background  to 

minimize  the  reflective  radiance  from  the  environment. 

1.  Amber  Engineering  AE4128 IR  Imaging  System 

The  Amber  system  is  a  128x128  Indium  Antimonide  staring  array,  which  operates 
in  the  3-5  jlm.  wavelength  band.  Its  display  provides  gray-scale  shading  or  pseudocolor. 

The  camera  assembly  has  a  cyrogenic  dewer,  which  is  filled  with  liquid  nitrogen  to  keep 
the  detector  temperature  around  77°K.  Its  capacity  is  400  ml.  To  provide  this  temperature 
level  and  to  prevent  the  condensation  on  the  detector  this  dewer  needs  to  be  evacuated  to 
less  than  lO'^  torr.  The  operator  can  control  contrast,  global  gain,  global  offset,  integration 
time,  and  frame  rate  to  adjust  the  camera  [Ref.  24].  In  this  experiment  global  gain  was  set 
to  one,  global  offset  was  adjusted  to  zero,  frame  rate  was  109  Hz,  and  integration  time 
was  62.  The  two  point  calibration  method  [Ref.  24]  was  used  to  calibrate  the  system. 
Figure  V.2  is  a  picture  of  the  Amber  camera. 

2.  Mitsubishi  Electronics  IR-M500  Thermal  Imager 

This  system  does  not  need  any  external  cooling  process.  The  Mitsubishi  camera 
has  a  cyrogenic  Stirling  cycle  cooler  to  keep  the  detector  temperature  around  80'’K.  This  is 
a  512x512  PtSi  staring  array,  which  operates  in  the  3-5  |lm.  wavelength  band  operator 
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has  a  small  wired  remote  controller  and  can  not  control  the  most  of  the  features.  This 
system' s  picture  is  shown  in  Figure  V.3. 


B.  PROCEDURE 


As  a  first  step  the  input  aperture  lens  of  the  TIS  and  the  target  bar  pattern  were 
adjusted  to  the  same  height  and  then  aligned  perpendicular  to  each  other.  As  a  second 
step  each  TIS  was  calibrated.  The  laboratory  was  darkened  during  the  measurements.  To 
familiarize  the  observers  with  four  bar  pattern' s  position  and  view  on  display  the 
background  plate  was  heated  up  and  observer  was  allowed  to  adjust  the  camera  and  his 
position  to  have  optimum  view.  During  the  measurements  observers  knew  where  to  look 
on  the  d  isplay  for  the  four  bar  pattern.  Measurements  were  done  in  two  stages  for  each 
system  with  multiple  observers.  In  the  first  stage  the  background  plate  was  heated  up 
until  the  observer  decided  that  four  bar  pattern  was  100%  resolvable,  then  the  temperature 
difference  between  the  thermocouple  readings  were  recorded  as  the  MRTD.  In  the  second 
stage,  the  backplate  was  cooled  down,  via  air  cooling,  until  the  observer  could  not  resolve 
100%  of  the  four  bar  pattern  and  again  the  temperature  difference  was  recorded  as  an 
MRTD.  These  steps  were  repeated  for  different  spatial  frequencies.  An  entire  set  of 
measurements  were  taken  by  one  observer ,  then  same  process  was  done  by  another 
observer.  This  process  was  repeated  for  both  horizontal  and  vertical  directions  in  order  to 
characterize  the  MRTD  values  in  these  directions  and  to  calculate  the  MRTDjd- 

C.  EXPERIMENTAL  MRTD  MEASUREMENTS 

MRTD  measurements  showed  that  MRTD^  has  a  higher  value,  or  equivalently  a 
lower  performance,  than  MRTD^  for  both  of  the  systems.  Figure  V.4.  shows  the  MRTDh, 
MRTD„  and  MRTD2d  curves  for  Amber  AE4128  IR  Imaging  System.  MRTDjd  was 
calculated  from  MRTDyand  MRTDh  using  Equation  (III.13).  Figure  V.5.  is  the 
logarithmic  presentation  of  the  same  values.  Figure  V.6  shows  the  MRTD^,  MRTDy,  and 
MRTD2d  values  for  Mitsubishi  IR-M500  Thermal  Imager.  Figure  V.7  presents  the  same 
values  in  logarithmic  scale  [Ref.  27]. 
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Figure  V.4.MRTD  Measurements  for  Amber  AE4128 IR  Imaging  System 
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Figure  V.5.  Logarithmic  Presentation  of  the  MRTD  Measurements  for 
Amber  AE4128  IR  Imaging  System 
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Spatial  Frequency  (cyc/mr) 

Figure  V.6.  MRTD  Measurements  for  Mitsubishi  IR-M500  Thermal  Imager 


D.  OBJECTIVE  MEASUREMENTS 


These  experiments  were  done  with  the  Mitsubishi  System  and  instead  of  a  human 
observer  an  oscilloscope  was  used  to  model  an  ATR.  Three  different  sets  of  measurement 
were  taken. 

1.  Constant  SNR 

For  objective  MRTD  measurements  the  acceptable  SNR  was  set  to  6.0  which 
gives  reasonable  target  resolution  on  the  measurement  device,  an  oscilloscope.  Using  the 
oscilloscope  AV signal  and  AV noise  were  measured  and  their  ratio  used  as  the  SNR.  When 
SNR  reached  6.0,  the  temperature  difference  between  the  front  and  back  plate  was 

recorded  for  four  bar  patterns  of  different  spatial  frequencies.  Figure  V.8  shows  that 
MRTD  measured  with  Tektronix  468  Digital  Storage  Oscilloscope.  This  value  for  MRTD 
is  much  higher  than  the  MRTD  measured  by  the  human  observer.  The  MRTD  value  can 
change  from  observer  to  observer,  but  the  objective  MRTD  is  more  reproducible. 


Figure  V.8.  Comparison  of  MRTD  Measured  by  ATR  to  MRTD  Measured 
by  Human  Observer  for  Mitsubishi  System 
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Figure  V.9.  presents  the  picture,  which  shows  the  four  bar  pattern  view  on  the 
display  and  ATR  unit  display,  when  spatial  frequency  is  0.8  cyc/mr,  temperature 
difference  between  the  plates  is  2.69  °C  and  SNR  is  6.0. 


a.  Target  View  on  the  Display  b.  ATR  Unit  Display 


Figure  V.9.  The  View  of  the  Target  on  the  Monitor  and  ATR  Unit  Display 

From  the  picture  of  the  display.  Figure  V.9b,  noise  level  is  measured  as  0.4  unit  and 
signal  level  is  measured  as  2.4  unit.  The  ratio  of  the  signal  level  to  the  noise  level  gives 
6.0  as  accepted.  The  defined  SNR  threshold  of  6.0  is  selected  arbitrarily. 

2.  Constant  Temperature 

In  a  second  set  of  measurements  the  temperature  difference  between  the  plates 
was  kept  constant,  1.94  °C,  and  SNR  was  measured  by  using  the  ATR  unit  for  different 
spatial  frequencies.  Figure  V.IO  presents  the  results  of  this  measurement.  Figure  V.8 
shows  that  a  1.92  °C  temperature  difference  between  the  plates  is  needed  to  have  SNR 
equal  to  6.0,  when  the  spatial  frequency  is  0.4  cyc/mr.  If  this  result  is  compared  to  the 
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Figure  V.  10,  it  is  seen,  as  expected,  that  SNR  starts  to  go  down  for  values  higher  than  0.4 
cyc/mr. 


Figure  V.IO  SNR  for  Different  Spatial  Frequencies  at  Constant 
Temperature,  1.92  °C,  by  Using  ATR 


3.  Constant  Spatial  Frequency 

In  a  third  set  of  measurements  the  spatial  frequency  was  chosen  constant  at  0.5 
cyc./mr.  and  a  temperature  difference  between  the  plates  was  recorded  for  the  various 
SNR  levels,  which  were  measured  by  ATR,  from  1.0  to  8.0.  The  results  of  the 
measurement  are  presented  in  Figure  V.l  1  As  expected  SNR  level  increases  temperature 
difference  between  the  plates  increases. 

These  three  objective  measurements  demonstrated  the  interesting  relationships 
between  the  SNR,  temperature  difference  between  the  plates,  and  spatial  frequency. 
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keeping  one  of  them  constant  each  time.  At  the  same  time  the  ATR  unit  MRTD  was 
compared  to  the  subjective  MRTD. 


SNR 

Figure  V.l  1 .  SNR  vs.  Delta  T  (°C)  at  0.5  cyc./mr  Spatial  Frequency 

This  comparison  (Figure  V.8)  shows  the  oscilloscope-measured  MRTD  to  exceed  the 
"subjective"  value  by  a  factor  of  approximately  5.  This  factor  is  strongly  influenced  by 
the  choice  of  "threshold"  SNR.  From  Figure  V.9.b  we  see  that  the  oscilloscope  clearly 
exceeds  the  noise:  this  must  be  considered  well  above  threshold.  Figure  V.9.a  is  also  seen 
to  be  clearly  resolved.  Thus  the  selection  6.0  clearly  exceeds  the  threshold  value  of  SNR 
needed  for  subjective  MRTD  evaluation,  and  the  objective  "MRTD"  curve  is  artificially 
elevated.  A  lower  threshold  SNR  selection  would  bring  the  "subjective"  and  "objective" 
curves  of  MRTD  into  closer  agreement.  However,  it  was  discovered  that  if  a  significantly 
lower  SNR  (i.e.,  SNR=2)  was  employed  in  the  oscilloscope  measurement  then  the  pattern 
appearing  on  the  oscilloscope  display  was  not  always  identifiable  as  a  four  bar  pattern. 
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VI.  COMPARISON  OF  FLIR92  MRTD  TO 
MEASURED  MRTD  AND  CONCLUSION 


A.  COMPARISON  OF  FLIR92  MRTD  TO  MEASURED  MRTD 

MRTD  measurements,  which  made  with  the  Amber  Engineering  AE4128  IR 
Imaging  System  and  Mitsubishi  Electronics  IR-M500  Thermal  Imager,  were  presented  in 
Chapter  V.  This  chapter  compares  the  measured  MRTD  values  to  the  predicted  FLIR92 
MRTD  values.  Short  output  listings  of  the  FLIR92  code  for  both  of  the  systems  are 
presented  in  Appendix  D. 

1.  Amber  Engineering  AE4128  IR  Imaging  System 

Predicted  FLIR92  MRTD^,  MRTD^,  and  MRTD^^  for  the  Amber  system  can  be 
seen  in  Figure  VI.  1.  Figure  VI.2  compares  the  FLIR92  and  measured  MRlDj^.  Figure 
VI.  3  is  the  logarithmic  presentation  of  the  same  comparison.  These  figures  show  that 
predicted  FLIR92  MRTD  is  more  optimistic  than  the  measured  MRTD  by  a  factor  of 
about  one  hundred.  Figure  VI.  1  shows  predicted  MRTD  of  the  functional  form  expected, 
and  similar  shape  to  the  measurements  shown  in  Figure  VI.2.  However  discrepancy  in 
magnitude  of  about  a  hundred  is  seen  between  prediction  and  measurement.  Reasons  for 
this  may  be  sought  in  the  input  to  prediction  and  in  the  experience  level  of  the  operators 
making  the  measurements.  The  FLIR92  MRTD  value  assumes  "experienced  operators", 
which  makes  it  unrealistic  for  our  measurements.  Additionally  the  Amber  system  is  an 
Indium  Antimonide  Focal  Plane  Array  sensor,  which  necessitates  input  of  spectral 
detectivity  for  InSb. 
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2.  Mitsubishi  Electronics  IR-M500  Thermal  Imager 


Predicted  FLIR92  MRTD^,  MRTD^,  and  MRTD2d  for  Mitsubishi  system  are 
presented  in  Figure  VI.4.  Figure  VI.5  is  the  comparison  of  the  FLIR92  MRTD2d  to  the 
measured  MRTD2d.  Figure  VL6  shows  this  comparison  in  logarithmic  scale.  Predicted 
FLIR92  MRTD  is  better  than  measured  MRTD,  but  these  values  are  closer  to  each  other 
than  Amber  system's  values.  The  corresponding  comparison  for  the  Mitsubishi  system 
(Figure  VI.5)  shows  much  closer  agreement,  with  again,  the  modeled  MRTD  being  more 
optimistic  (lower)  than  the  measured.  The  discrepancy  is  now  in  the  range  of  50%  of 
predicted,  which  might  be  accounted  for  by  the  experience  factor  in  the  operators. 


0  0.2  0.4  0.6  0.8  1 

Spatial  Frequency  (cyc./mr) 

Figure  VI.  1.  FLIR92  MRTD  Values  for  Amber  System 
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Spatial  Frequency  (cyc./mr) 


Figure  VI.2.  Comparison  of  the  FLIR92  MRTD  to  the  Measured  MRTD 

for  Amber  System 


Spatial  Frequency  (cyc./mr) 

Figure  VI.  3.  Logarithmic  Comparison  of  the  FLIR92  MRTD  to  the 
Measured  MRTD  for  Amber  System 


Spatial  Frequency  (cyc./mr) 


Figure.  VI.4.  FLIR92  MRTD  Values  for  Mitsubishi  System 


Spatial  Frequency  (cyc./mr) 

Figure  VI.5.  Comparison  of  the  FLIR92  MRTD  to  the  Measured  MRTD 

for  Mitsubishi  System 


Figure  VI.6.  Logarithmic  Comparison  of  the  FLIR92  MRTD  to  the 
Measured  MRTD  for  Mitsubishi  System 


B.  LIMITATIONS  OF  THE  FLIR92  AND  CONCLUSION 

1.  FLIR92  TIS  Performance  Model's  3D  noise  concept  was  covered  in  Chapter  TTT 
Only  NVEOD  has  the  necessary  technology  to  measure  these  3D  noise  components.  For 
this  reason  default  noise  values  were  used  in  the  applications. 

2.  FLIR92  does  not  cover  the  aliasing  effect.  For  this  reason  it  does  not  make  any 
MRTD  prediction  for  spatial  frequencies  higher  than  the  Nyquist  Frequency.  The  aliasing 
effect  was  mentioned  in  Chapter  IV,  FLIR92  needs  to  include  this  effect  to  make  MRTD 
predictions  for  spatial  frequencies  higher  than  the  Nyquist  Frequency. 

3.  FLIR92  does  not  predict  the  MRTD  correctly  at  low  spatial  frequencies,  its 
value  is  optimistic  [Ref.  2  .p.  397]. 

4.  The  FLIR92  validation  data  base  includes  only  PtSi  and  HgCdTe  sensors  for 
staring  systems.  For  this  reason  the  detectivity  curve  for  Amber  Indium  Antimonide 
staring  array  was  identified  by  using  spectral  detectivity  data  at  three  wavelength  points. 
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The  Mitsubishi  PtSi  staring  array  system  did  not  need  to  be  defined  by  spectral  detectivity 
data.  Comparisons  of  the  FLIR92  MRTD  values  to  the  measured  MRTD  values  show 
that  there  is  big  difference  between  measured  and  predicted  MRTD  for  the  Amber 
System,  whereas  the  values  for  the  Mitsubishi  System  are  close  together.  If  Indium 
Antimonide  is  added  to  the  FLIR92  validation  data  base,  the  MRTD  predictions  for  this 
system  may  be  closer  to  the  measured  values. 

5.  Explanations  and  equations,  which  are  given  in  Ref.(9)  and  (10),  are  not 
adequate  to  describe  me  FLIR92  TIS  Performance  Model  exacUy.  Source  code  for  this 

T',odel  is  needed  to  make  this  system  concept  clearer  and  to  prepare  a  background  for 
future  studies. 

6.  FLIR92  does  not  plot  the  predicted  output  values.  A  plot  program  is  needed  to 
plot  the  outputs. 
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APPENDIX  A.  THERMAL  IMAGING  SYSTEM  PARAMETERS 


A.  SAMPLE  FIRST  GENERATION  SYSTEM 

In  Chapter  n  to  present  Ratches  Model,  Lloyd  Model,  and  FLIR92  Model  this 
sample  system  was  chosen  in  similarity  with  the  one  employed  by  Lloyd  [Ref.  1]  and 

used  in  Reference  3.  The  parameters  of  this  system  are  as  follows: 

Lens  focal  length  (d) . 50000pm 

Diameter  of  lens  (D) . 20000pm 

F/  number . 2.5 

Detector  array  individual  element  size  (square) . 0.005cm 

oc . Imrad 

P . Imrad 

Detectors  cold  shielding . not  background  limited 

Characteristic  wavelength  of  the  detectors . 1 1.5pm 

Spectral  bandpass  of  the  detectors . 8pm-  1 1.5pm 

Specific  detectivity . 2x10'°  cm  Hz°  Vwatt 

Frame  rate . 30  Hz 

Scan  rate  format . 60  fields/sec 

Number  of  detectors  in  parallel . 150 

Number  of  scan  lines . 300 

Interlace . 2  to  1 

Horizontal  scan  efficiency . 0.8 

Vertical  scan  efficiency . 0.8 

Overall  scan  efficiency . 0.64 

Overscan  ratio . 1 

Distance  between  horizontal  scan  lines . . . Imrad 

HFOV . 400mrad 

VFOV.... . SOOmrad 

Detector  dwell  time . 2.67x10'^  sec 
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Horizontal  scanning  velocity . 

3-dB  freq.  electronic  roll-off. . 

Noise  equivalent  reference  bandwidth . 

SNRj^  for  detection  of  one  bar . 

Background  temperature . 

Monochromatic  wavelength  of  the  target. 
Optical  efficiency  of  the  viewer . 


.37453.2mrad/sec 

. 18716.6  Hz 

. 29.4KHZ 

. 4.5 

. 300°K 

. lOpm 

. 0.8 


B.  SECOND  GENERATION  SYSTEMS 


The  Amber  Engineering  AE4128  IR  Imaging  System  and  Mitsubishi  Electronics 
IR-M500  Thermal  Imager  were  used  for  laboratory  measurements  and  their  results  were 
compared  to  predicted  FLIR92  values.  The  parameters  of  these  systems  are  as  follows: 

1.  Amber  Engineering  AE4128  IR  Imaging  System 


BLIP  Performance . 

Spectral  cut-on . 

Spectral  cut-off. . 

F/  number . 

Focal  length . 

Optical  transmittance . 

Frame  rate . 

Detector  active  horizontal  dimension 
Detector  active  vertical  dimension 

D* . : . ” 

Integration  time . 

Number  of  horizontal  detectors . 

Number  of  vertical  detectors . 

Detector  cell  horizontal  dimension . 


. YES 

. S.Ojim 

. 5.0  |im 

. 3.0 

. 10.0  cm 

. 0.95 

. 109  Hz 

. 40  |im 

. 40)i.m 

5.97x10”  cm  Hz°  Vwatt 

. 8887.615  |i.sec 

. 128 

. 128 

. 50  lim 
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Detector  cell  vertical  dimension 
Number  of  active  CRT  lines . 

Display  brightness . 

3D  noise  level . 


. 50  [im 

. 480 

10.0  mLamberts 
. MOD 


2.  Mitsubishi  Electronics  IR-M500  Thermal  Imager 


BLIP  performance . YES 

Spectral  cut-on . 3.0p.m 

Spectral  cut-off. . 5.0|im 

F/ number . 1.4 

Focal  length . 5.0  cm 

Optical  transmittance . 0.95 

Frame  rate . 60  Hz 

Detector  active  horizontal  dimension . 16.24  p.m 

Detector  active  vertical  dimension . 12.49}im 

D* . . . 5.0xl0’°cm  Hz°  Vwatt 

Integration  time . 16145.833 

Number  of  horizontal  detectors . 512 

Number  of  vertical  detectors . 512 

Detector  cell  horizontal  dimension . 26  (im 

Detector  cell  vertical  dimension . 20  (im 

Number  of  active  CRT  lines . 480 

PtSi  emission  coefficient . 0.16  1/eV 

Schottky  barrier  height . 22  eV 

Display  brightness . 10.0  mLamberts 

3D  noise  level . MOD 
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APPENDIX  B.  FLIR92  MTF  EQUATIONS  AND  D  OPERATOR 


FLIR92  TIS  Performance  Model  was  covered  in  Chapter  III.  This  appendix  gives 
all  MTF  equations  used  in  FLIR92  ( Ref.  [9] )  and  a  different  approach  to  the  usage  of  D 
operators. 


A.  MTF  EQUATIONS 

1.  Prefilter  MTFs 

a.  Optics  MTFs 

(1)  Diffraction-limited  MTF 


(2)  Geometric  Blur  MTF 


b.  Detector  Spatial  MTF 


sin(ji5zf) 

7c8zf 


(A.3) 


c.  Focal  Plane  Array  Integration  Time 


H  (f)  ^  linWvxti) 


(A.4) 
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d.  Sample-scene  Phase  MTF 


Hssp(f)  =  COS 


(A.5) 


e.  Image  Motion  MTFs 

(1)  Linear  Image  Motion  MTF 


sin(7tVrtif) 
71  Vr  tjf 


(A.6) 


(2)  Random  Image  Motion  MTF 


Hmr(f)  =  e 


-2k^  r 


(A.7) 


(3)  Sinusoidal  Image  Motion  MTF 


H„(f)  =  Jo(2jtAf) 


(A.8) 


2.  Temporal  Postfilter  MTFs 

a.  Detector  Temporal  MTF 


Hdt(ft)=  1  + 


(A.9) 


b.  Electronics  Low  Frequency  Response 


^ehp 


(A.  10) 
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c.  Electronics  High  Frequency  Response 


Help(ft)  = 


f 

1  + 


2n  \-0-5 


(A.  11) 


d.  Boosting  MTF 


(A.  12) 


3.  Spatial  Postfilter  MTFs 

a.  Electro-optical  Multiplexor  MTF 


Heom(f)  = 


sin(7tf5ied) 


7Cf5 


led 


(A.  13) 


b.  Digital  Filter  MTF 

H<„g(f)=  S  aicosj^^j 

Hd,g(f)  =  2  aicos 

c.  Display  MTF 

(1)  CRT  Display  MTF 

( 

-log(0.025) 


if  N  is  odd 


if  N  is  even 


(A.14) 


(A.  15) 


(A.  16) 


(A.17) 
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d.  CCD  Charge  Transfer  Efficiency  MTF 


e.  Display  Sample  and  Hold  MTF 


Hdsh(f) 


sin(i:5sf) 

Jt5sf 


(A.18) 


(A.19) 


f.  Eye  MTF 

(1)  Non-limiting  Eye  MTF 


Heye(f)=1.0 


(A.20) 


(2)  Limiting  Eye  MTF 


(A.21) 


B.  DIRECTIONAL  AVERAGING  OPERATOR  (  D  OPERATOR) 


Reference  12  brings  another  approach  to  the  usage  of  D  operators.  In  this 
approach  by  applying  D,,  D^,  or  D„  operations  to  the  noise  components  and  by  using 
subtraction  between  the  noise  data  sets,  we  extract  a  desired  noise  component.  As  an 
example,  if  we  want  to  find  N^, ,  the  following  steps  must  be  applied,  which  are  shown  in 
Figure  III.  3; 

a.  Apply  Dj  and  D^  to  the  composite  data  set: 

(1)  D(  process  gives  data  set-1:  ,  N^,  N^,  and  s 

(2)  D^  process  gives  data  set-2:  N^,  N^,  N^,  and  s 

b.  Apply  D^  to  the  data  set-1  to  obtain  the  data  set-3:  and  s 
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c.  Subtract  the  data  set-3  from  the  data  set-2  to  have  data  set-4:  and  s 

d.  Apply  D(  to  the  data  set-3  to  obtain  the  s-component 

e.  Subtract  s-component  ( step,  d )  from  data  set-4.  This  gives  the  final 
data  set  which  has  only  Njj, . 


Figure  B.l.  Application  of  D  Operator 
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APPENDIX  C.  LABORATORY  MEASUREMENTS 


This  appendix  presents  the  measured  data,  which  is  mentioned  in  Chapter  V  and 
shows  plots  of  MRTD  measured  during  the  heating  and  cooling  cycles  in  horizontal  and 
vertical  directions  for  the  Amber  and  Mitsubishi  systems. 


Table  C.l.  Amber  System  Measured  MRTD,j 


Spatial  Freq. 
(cyc./mr) 

MRTD-Heating 

(°C) 

MRTD-Cooling 

(°C) 

Average  MRTD 
CO 

0.1 

0.33 

0.28 

0.31 

0.15 

0.33 

0.31 

0.32 

0.2 

0.33 

0.17 

0.25 

0.25 

0.28 

0.39 

0.34 

0.3 

0.36 

0.31 

0.34 

0.375 

0.33 

0.44 

0.39 

0.4 

0.39 

0.28 

0.34 

0.50 

0.33 

0.28 

0.31 

0.60 

0.44 

0.44 

0.44 

0.625 

0.42 

0.33 

0.38 

0.70 

0.56 

0.56 

0.56 

0.75 

0.61 

0.61 

0.61 

0.80 

0.58 

0.75 

0.67 

0.875 

0.67 

0.83 

0.75 

1.0 

0.78 

1.6 

1.11 

1.2 

1.5 

2.6 

2.05 
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Table  C.2.  Amber  System  Measured  MRTD^ 


Spatial  Freq. 
(cyc./mr) 

MRTD-Heating 

(°C) 

MRTD-Cooling 

(°C) 

Average  MRTD 
CO 

0  10 

0.33 

0.28 

0.31 

0.15 

0.33 

— 

0.31 

0.32 

!  0.20 

0.33 

0.17 

0.25 

i  0.25 

0.28 

0.39 

0.33 

0.30 

0.36 

0.31 

0.33 

0.375 

0.39 

0.39 

0.39 

0.40 

0.44 

0.44 

0.44 

0.50 

0.19 

0.56 

0.38 

0.60 

0.25 

0.64 

0.44 

0.625 

0.42 

0.61 

0.51 

0.70 

0.56 

0.81 

0.68 

0.75 

1.1 

0.89 

1.0 

0.80 

1.1 

0.83 

0.94 

0.875 

1.6 

1.0 

1.3 

1.0 

1.6 

I.l 

1.4 

1.2 

3.3 

3.3 

3.3 
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Tabl 

e  C.3.  Mitsubishi 

System  Measured  Ml 

Spatial  Freq. 
(cyc./mr) 

MRTD-Heating 

(“C) 

MRTD-Cooling 

ro 

Average  MRTD 
CO 

0.10 

0.08 

0.06 

0.07 

0.15 

0.08 

0.06 

0.07 

0.20 

0.12 

0.18 

0.15 

0.25 

0.15 

0.13 

0.14 

0.30 

0.28 

0.11 

0.19 

0.375 

0.26 

0.13 

0.19 

0.40 

0.28 

0.14 

0.21 

0.50 

0.26 

0.13 

0.19 

0.60 

0.28 

0.19 

0.24 

0.625 

0.31 

0.19 

0.25 

0.70 

0.36 

0.31 

0.33 

0.75 

0.28 

0.25 

0.26 

0.80 

0.36 

0.36 

0.36 

0.875 

0.31 

0.31 

0.31 

1.0 

0.44 

0.33 

0.39 

1.2 

0.78 

0.78 

0.78 
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Tab 

le  C.4.  Mitsubishi  System  Measured  MRTD„ 

Spatial  Freq. 
(cyc./mr) 

MRTD-Heating 

CO 

MRTD-Cooling 

(“C) 

Average  MRTD 
CC) 

0.10 

0.08 

0.06 

0.07 

0.15 

0.08 

0.06 

0.07 

0.20 

0.17 

0.11 

0.14 

0.25 

0.22 

0.11 

0.17 

0.30 

0.25 

0.11 

0.18 

0.375 

0.22 

0.14 

0.18 

0.40 

0.39 

0.11 

0.25 

0.50 

0.39 

0.19 

0.29 

0.60 

0.44 

0.31 

0.38 

0.625 

0.44 

0.47 

0.46 

0.70 

0.81 

0.44 

0.63 

0.75 

0.83 

0.75 

0.79 

0.80 

1.1 

0.70 

0.90 

0.875 

1.3 

0.69 

1.3 

1.0 

* 

* 

1.2 

* 

* 

* 
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Table  C.5.  Constant  SNR=  6.0 


Spatial  Freq. 
(cyc./mr) 

Delta  -T 
(°C) 

0.2 

1.61 

0.3 

1.61 

0.4 

1.86 

0.5 

2.0 

0.6 

2.36 

0.7 

2.39 

0.8 

2.69 

Table  C.6.  Constant  Delta-T  =1.94°C 


Spatial  Freq. 
(cyc./mr) 

SNR 

0.2 

6.25 

0.4 

6.23 

0.5 

5.5 

0.6 

4.75 

0.7 

4.25 

0.8 

3.38 

Table  C.7.  Constant  Spatial  Freq.=  0.5  cyc./mr 


SNR 

Delta-T 

rc) 

1 

0.8 

2 

1.5 

3 

1.8 

4 

2.35 

5 

3.0 

6 

3.7 

7 

4.2 

8 

5.0 
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0 


0.5 


1 


Spatial  Frequency  (cvc/mr) 


Figure  C.l.  Amber  System  Measured  MRTD^ 


Spatial  Frequency  (cyc/mr) 


Figure  C.2.  Amber  System  Measured  MRTD^ 
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Spatial  Frequency  (cyc/mr) 


Figure  C.3.  Mitsubishi  System  Measured  MRTD^ 


Spatial  Frequency  (cyc/mr) 

Figure  C.4.  Mitsubishi  System  Measured  MRTD^ 
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APPENDIX  D.  FLIR92  SHORT-LISTING  OUTPUTS 


This  Appendix  shows  the  short-listing  outputs  of  the  FLIR92  for  Amber 
Engineering  AE4128  IR  Imaging  System  and  Mitsubishi  Electronics  IR-M500  Thermal 
Imager. 


A.  AMBER  ENGINEERING  AE4128  IR  IMAGING  SYSTEM 


U.S.  Army  CECOM  NVESD  FLIR92  Wed  Sep  27  20:01:16  1995 

output  file:labout  1 . 1  short  listing 
data  file:  labinl 

command  line  arguments;  -d  labinl  -o  laboutl  -p  BOTH  -a  laboutl 

begin  data  file  listing  .  . . 
flir92  data  file  template 
>environment 


laboratory_temperature 

300.0 

K 

background_temperature_  1 

294.11 

K 

BLIP_performance 

YES 

YES_or 

>spectral 

spectral_cut_on 

3.0 

microns 

spectral_cut_off 

5.0 

microns 

diffraction_wavelength 

0.0 

microns 

>optics_l 

eff_f_number 

3.0 

— 

eff_focal_length 

10.0 

cm 

efif_aperture_diameter 

0.0 

cm 

optics_blur_spot 

0.0 

mrad 

average_optical_trans 

0.0 

>optics_2 

HFOV:  VFOV_aspect_ratio 

0.0 

— 

magnification 

0.0 

— 

frame_rate 

109.0 

Hz 

fields_per_frame 

0.0 

— 
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xietector 


horz_dimension_(active) 
vert_dimension_(active) 
peak;_D_star 
integration_time 
l/f_knee_frequency 
>^a_stare 
#_horz_detectors 
#_vert_detectors 
horz_unit_cell_dimension 
vert_unit_ceU_dimension 
>scene_phasing 
horz_target/detector_phase 
vert_target/detector_phase 
>crt_display 
#_active_lines_on_CRT 
horz_cit_spot_sigma 
vert_crt_spot_sigma 
>electronics 
high_pass_3db_cuto  n 
high_pass_filter_order 
low_pass_3db_cutoff 
low_pass_filter_order 
boost_amplitude 
boost_frequency 
sample_and_hold 
>display 

display_brightness 

display_height 

display_viewing_distance 

>eye 

threshold_SNR 

eye_integration_time 

MTF 

>3d_noise_default 

noise_level 

>spectral_detectivity 


40.0  microns 
40.0  microns 
5.97ell  cm-sqrt(Hz)AV 
8887.615  microsec 
0.0  Hz 

128.0  - 
128.0 

50.0  microns 
50.0  microns 

0.0  degrees 
0.0  degrees 

480.0 

0.0  mrad 
0.0  mrad 

0.0  Hz 
0.0 

0.0  Hz 
0.0 
0.0 

0.0  Hz 

HORZ  HORZ_VERT_or_NO 

10.0  miBi-Lamberts 
17.78  cm 
88.9  cm 

6.0 

0.1  sec 

EXP  EXP_or_NL 
LO  NO_LO_MOD_or_HI 


#_points:  3  microns _ detectivity 

3.0  0.67 

3.34  0.8375 

4.5  1.0 
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>end 

end  data  file  listing  . . . 

MESSAGES 

diagnosticO:  Using  default  3D  noise  components. 
diagnostic():  Using  _LO_  level  3D  noise  defaults. 
diagnosticO:  Diffraction  wavelength  set  to  spectral  band  midpoint. 
diagnosticO:  Optics  transmittance  defaulted  to  0.7. 
diagnosticO:  Interlace  (fields  per  frame)  defaulted  to  2. 
diagnosticO:  Fields-of-view  calculated  by  model. 


CALCULATED  SYSTEM  PARAMETERS 

field-of-view:  3.666h  x  3.666v  degrees 
63.98h  X  63.98v  mrad 
magnification:  3.116 

optics  blur  spot:  29.280  microns  (diffraction-bmited) 

0.293  mrad 


detector  IFOV: 
FPA  fill  factor: 
FPA  duty  cycle: 


0.400h  X  0.400v  mrad 
0.640 
1.938 


TEMPERATURE  DEPENDENCE 
BLIP  detector 
scaling  factors  (T 1:300) 

NETD:  1.07 

peakD-star:  1.11 

Planck  thermal  derivative:  0.83 


parameter  NETD  @  300  K 


white  NETD 
classical  NETD 
sigma_TVH  NETD 
sigma.VH  NETD 


NETD  @  294  K  noise 
bandwidth 


0.010  deg  C 
0.010  deg  C 
0.008  deg  C 
0.003  deg  C 


0.010  deg  C 
0.010  deg  C 
0.008  deg  C 
0.003  deg  C 


8.561e4001  Hz 
8.565e+001  Hz 
5.629e+001  Hz 
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Planck  integral  2.127e-005  1.776e-005  W/(cm*cm*K) 

.  ..w/D-star  1.251e+007  1.166e-K)07  sqrt(Hz)/(cm*K) 


TOTAL  HORIZONTAL  MTFs 


cy/mr 

H_SYS 

H_PRE 

H_TPF 

H_SPF 

0.000 

1.000 

1.000 

1.000 

1.000 

0.125 

0.928 

0.977 

1.000 

0.950 

0.250 

0.841 

0.946 

1.000 

0.889 

0.375 

0.744 

0.908 

1.000 

0.819 

0.500 

0.642 

0.864 

1.000 

0.743 

0.625 

0.540 

0.814 

1.000 

0.663 

0.750 

0.441 

0.760 

1.000 

0.581 

0.875 

0.350 

0.702 

1.000 

0.499 

1.000 

0.269 

0.641 

1.000 

0.420 

1.125 

0.199 

0.579 

1.000 

0.344 

1.250 

0.141 

0.515 

1.000 

0.274 

1.375 

0.095 

0.452 

1.000 

0.210 

1.500 

0.060 

0.390 

1.000 

0.153 

1.625 

0.034 

0.329 

1.000 

0.103 

1.750 

0.017 

0.270 

1.000 

0.062 

1.875 

0.006 

0.215 

1.000 

0.027 

2.000 

0.000 

0.163 

1.000 

0.000 

2.125 

0.000 

0.115 

1.000 

-0.021 

2.250 

0.000 

0.072 

1.000 

-0.035 

2.375 

0.000 

0.034 

1.000 

-0.045 

2.500 

0.000 

0.000 

1.000 

0.000 

TOTAL  VERTICAL  MTFs 

cy/mr 

H_SYS 

H_PRE 

H_SPF 

0.000 

1.000 

1.000 

1.000 

0.125 

0.934 

0.977 

0.956 

0.250 

0.863 

0.946 

0.912 

0.375 

0.789 

0.908 

0.868 

0.500 

0.713 

0.864 

0.825 

0.625 

0.637 

0.814 

0.783 

0.750 

0.563 

0.760 

0.740 

0.875 

0.491 

0.702 

0.699 

1.000 

0.423 

0.641 

0.659 

1.125 

0.359 

0.579 

0.620 
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1.250 

0.300 

0.515 

0.582 

1.375 

0.246 

0.452 

0.545 

1.500 

0.198 

0.390 

0.509 

1.625 

0.156 

0.329 

0.475 

1.750 

0.120 

0.270 

0.442 

1.875 

0.088 

0.215 

0.411 

2.000 

0.062 

0.163 

0.381 

2.125 

0.041 

0.115 

0.352 

2.250 

0.023 

0.072 

0.325 

2.375 

0.010 

0.034 

0.300 

2.500 

0.000 

0.000 

0.276 

PREHLTER 

VALUES  AT  NYQUIST 

horz  H_PRE( 1 .00)  =  0.64 1  vert  H_PRE( 1 .00)  =  0.64 1 


SAMPLING  RATES 

horizontal  2.00  samples/mr 

vertical  2.00  samples/mr 

effective  2.00  samples/mr 

SENSOR  LIMITING  FREQUENCIES 


spatial  Nyquist 

horizontal  2.50  1.00 

vertical  2.50  1.00 

effective  2.50  1.00 

MRTD  3D  NOISE  CORRECTION  (AVERAGE) 

300  K  294  K 

horizontal  1.657  1.585 

vertical  1.657  1.585 


MRTD  AT  300  K  BACKGROUND  TEMPERATURE 


cy/mr 

horz 

cy/mr 

vert 

cy/mr 

2D 

0.05 

0.125 

0.001 

0.05 

0.125 

0.001 

0.127 

0.001 

0.10 

0.250 

0.003 

0.10 

0.250 

0.003 

0.150 

0.002 

0.15 

0.375 

0.005 

0.15 

0.375 

0.004 

0.174 

0.002 

0.20 

0.500 

0.007 

0.20 

0.500 

0.006 

0.197 

0.002 

0.25 

0.625 

0.010 

0.25 

0.625 

0.008 

0.221 

0.002 

0.30 

0.750 

0.013' 

0.30 

0.750 

0.011 

0.244 

0.003 

0.35 

0.875 

0.019 

0.35 

0.875 

0.014 

0.277 

0.003 

0.40 

1.000 

0.027 

0.40 

1.000 

0.018 

0.314 

0.004 
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0.45 

1.125 

99.999 

0.45 

1.125 

99.999 

0.351 

0.004 

0.50 

1.250 

99.999 

0.50 

1.250 

99.999 

0.391 

0.005 

0.55 

1.375 

99.999 

0.55 

1.375 

99.999 

0.439 

0.005 

0.60 

1.500 

99.999 

0.60 

1.500 

99.999 

0.486 

0.006 

0.65 

1.625 

99.999 

0.65 

1.625 

99.999 

0.538 

0.007 

0.70 

1.750 

99.999 

0.70 

1.750 

99.999 

0.591 

0.008 

0.75 

1.875 

99.999 

0.75 

1.875 

99.999 

0.646 

0.009 

0.80 

2.000 

99.999 

0.80 

2.000 

99.999 

0.702 

0.011 

0.85 

2.125 

99.999 

0.85 

2.125 

99.999 

0.762 

0.012 

0.90 

2.250 

99.999 

0.90 

2.250 

99.999 

0.823 

0.014 

0.95 

2.375 

99.999 

0.95 

2.375 

99.999 

0.882 

0.016 

1.00 

2.500 

99.999 

1.00 

2.500 

99.999 

0.936 

0.019 

MRTD 

AT  294  K  BACKGROUND 

TEMPERATURE 

cy/mr 

horz 

cy/mr 

vert 

cy/mr 

2D 

0.05 

0.125 

0.001 

0.05 

0.125 

0.001 

0.127 

0.001 

0.10 

0.250 

0.003 

0.10 

0.250 

0.003 

0.150 

0.002 

0.15 

0.375 

0.005 

0.15 

0.375 

0.005 

0.174 

0.002 

0.20 

0.500 

0.007 

0.20 

0.500 

0.006 

0.197 

0.002 

0.25 

0.625 

0.010 

0.25 

0.625 

0.009 

0.221 

0.002 

0.30 

0.750 

0.014 

0.30 

0.750 

0.011 

0.244 

0.003 

0.35 

0.875 

0.019 

0.35 

0.875 

0.015 

0.277 

0.003 

0.40 

1.000 

0.028 

0.40 

1.000 

0.019 

0.314 

0.004 

0.45 

1.125 

99.999 

0.45 

1.125 

99.999 

0.351 

0.004 

0.50 

1.250 

99.999 

0.50 

1.250 

99.999 

0.391 

0.005 

0.55 

1.375 

99.999 

0.55 

1.375 

99.999 

0.439 

0.006 

0.60 

1.500 

99.999 

0.60 

1.500 

99.999 

0.486 

0.006 

0.65 

1.625 

99.999 

0.65 

1.625 

99.999 

0.538 

0.007 

0.70 

1.750 

99.999 

0.70 

1.750 

99.999 

0.591 

0.008 

0.75 

1.875 

99.999 

0.75 

1.875 

99.999 

0.646 

0.010 

0.80 

2.000 

99.999 

0.80 

2.000 

99.999 

0.702 

0.011 

0.85 

2.125 

99.999 

0.85 

2.125 

99.999 

0.762 

0.013 

0.90 

2.250 

99.999 

0.90 

2.250 

99.999 

0.823 

0.015 

0.95 

2.375 

99.999 

0.95 

2.375 

99.999 

0.882 

0.017 

1.00 

2.500 

99.999 

1.00 

2.500 

99.999 

0.936 

0.019 

MDTD  AT  300  K  BACKGROUND  TEMPERATURE 


1/mr 

MDTD 

0.20 

12.500 

1.232 

0.40 

6.250 

0.315 

0.60 

4.167 

0.145 

0.80 

3.125 

0.086 
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1.00 

2.500 

0.058 

1.20 

2.083 

0.043 

1.40 

1.786 

0.034 

1.60 

1.563 

0.028 

1.80 

1.389 

0.023 

2.00 

1.250 

0.020 

2.20 

1.136 

0.018 

2.40 

1.042 

0.016 

2.60 

0.962 

0.014 

2.80 

0.893 

0.013 

3.00 

0.833 

0.012 

3.20 

0.781 

0.011 

3.40 

0.735 

0.010 

3.60 

0.694 

0.010 

3.80 

0.658 

0.009 

4.00 

0.625 

0.009 

4.20 

0.595 

0.008 

4.40 

0.568 

0.008 

4.60 

0.543 

0.007 

4.80 

0.521 

0.007 

5.00 

0.500 

0.007 

MDTD  AT  294  K  BACKGROUND  TEMPERATURE 
l/mr  MDTD 


0.20 

12.500 

1.266 

0.40 

6.250 

0.324 

0.60 

4.167 

0.149 

0.80 

3.125 

0.088 

1.00 

2.500 

0.060 

1.20 

2.083 

0.044 

1.40 

1.786 

0.035 

1.60 

1.563 

0.028 

1.80 

1.389 

0.024 

2.00 

1.250 

0.021 

2.20 

1.136 

0.018 

2.40 

1.042 

0.016 

2.60 

0.962 

0.015 

2.80 

0.893 

0.013 

3.00 

0.833 

0.012 

3.20 

0.781 

0.011 

3.40 

0.735 

0.011 
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3.60 

0.694 

0.010 

3.80 

0.658 

0.009 

4.00 

0.625 

0.009 

4.20 

0.595 

0.008 

4.40 

0.568 

0.008 

4.60 

0.543 

0.008 

4.80 

0.521 

0.007 

5.00 

0.500 

0.007 

FLIR92.  . .  laboutl.l;  end  of  listing 


B.  MITSUBISHI  ELECTRONICS  IR-MSOO  THERMAL  IMAGER 


U.S.  Army  CECOM  NVESD  FLIR92  Thu  Sep  28  09: 1 1  ;44  1995 

output  file:  mitsh.  1  short  listing 
data  file:  mitsh 

command  line  arguments:  -d  mitsh  -o  mitsh  -p  BOTH  -a  mitsh 

begin  data  file  listing  . . . 
flir92  data  file  template 
>environment 


laboratory_temperature 

300.0 

K 

background_temperature_  1 

294.0 

K 

BLIP_performance 

YES 

YES  or 

>spectral 

spectral_cut_on 

3.0 

microns 

spectral_cut_off 

5.0 

microns 

di£fraction_wavelength 

0.0 

microns 

>optics_l 

efif_f_number 

1.4 

eff_focal_length 

5.0 

cm 

eff_aperture_diameter 

0.0 

cm 

optics_blur_spot 

0.0 

mrad 

average_optical_trans 

0.95 

— 

>optics_2 

HFOV :  VFOV_aspect_ratio 

0.0 

_ 

magnification 

0.0 

_ 
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frame_rate 

fields_per_fraine 

xietector 

horz_diraension_(  active) 
vert_dimension_(active) 
peak_D_star 
integration_time 
l/f_knee_frequency 
>^a_stare 
#_horz_detectors 
#_vert_detectors 
horz_unit_cell_ditnension 
vert_unit_cell_dimension 
>PtSi 

emission_coefficient 

Schottky_banier_height 

>electronics 

high_pass_3db_cuton 

high_pass_filter_order 

low_pass_3db_cutoff 

low_pass_filter_order 

boost_amplitude 

boost_frequency 

sample_and_hold 

>display 

display_brightness 

display_height 

display_viewing_distance 

>crt_display 

#_active_lines_on_CRT 

horz_crt_spot_sigma 

vert_crt_spot_sigma 

>eye 

threshold_SNR 

eye_integration_time 

MTF 

>3d_noise_default 

noise_level 

>end 

end  data  file  listing  . . . 


60.0  Hz 

1.0  - 

16.24  microns 
12.49  microns 
5.0el0  cm-sqrt(Hz)AV 
16145.833  microsec 
0.0  Hz 

512  - 

512 

26  microns 
20  microns 

0.16  1/eV 
0.22  eV 

0.0  Hz 
0.0  - 
0.0  Hz 
0.0  - 
0.0  - 
0.0  Hz 

HORZ  HORZ_VERT_or_NO 

10.0  milli-Lamberts 
27.94  cm 
88.9  cm 

480 

0.0  mrad 
0.0  mrad 

6.0  - 

0.2  sec 

EXP  EXP_or_NL 

MOD  NO_LO_MOD_or_HI 
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MESSAGES 

diagnostic():  Using  default  3D  noise  components. 
diagnosticQ:  Using  _LO_  level  3D  noise  defaults. 
diagnosticO:  Diffraction  wavelength  set  to  spectral  band  midpoint. 
diagnostic();  Fields-of-view  calculated  by  model. 
diagnosticO:  PtSi  spectral  detectivity  predicted  by  model. 


CALCULATED  SYSTEM  PARAMETERS 

field-of-view:  15. 166h  x  1  L694v  degrees 
264.68h  X  204.09v  mrad 
magnification:  1.527 

optics  blur  spot:  13.664  microns  (diSraction-hmited) 

0.273  mrad 

detector  IFOV:  0.325h  x  0.250v  mrad 
FPA  fill  factor:  0,390 

FPA  duty  cycle:  0.969 


NORMALIZED  DETECTOR  SPECTRAL  DETECTIVITY 
wavelength  detectivity 


3.00 

1.00 

3.22 

0.78 

3.44 

0.60 

3.67 

0.46 

3.89 

0.34 

4.11 

0.24 

4.33 

0.17 

4.56 

0.11 

4.78 

0.07 

5.00 

0.03 

TEMPERATURE  DEPENDENCE 
BLIP  detector 
scaling  factors  (T 1:300) 

NETD:  1.10 

peakD-star:  1.12 

Planck  thermal  derivative:  0.83 
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parameter  NETD  @  300  K 


NETD  @  294  K 


noise 

bandwidth 


white  NETD  0.200  deg  C 

classical  NETD  0.200  deg  C 
sigma_TVH  NETD  0. 162  deg  C 


0.221  degC 
0.221  deg  C 
0.179  deg  C 


sigma_VH  NETD  0.065  deg  C 


0.065  deg  C 


Planck  integral  2. 1 27e-005 

. . .  w/D-star  1.986e-i-005 


1.769e-005 

1.802e+005 


TOTAL  HORIZONTAL  MTFs 


cy/rnr 

H_SYS 

H_PRE 

H_TPF 

H_SPF 

0.000 

1.000 

1.000 

1.000 

1.000 

0.154 

0.850 

0.974 

1.000 

0.872 

0.308 

0.679 

0.940 

1.000 

0.721 

0.462 

0.509 

0.900 

1.000 

0.566 

0.616 

0.358 

0.853 

1.000 

0.419 

0.770 

0.235 

0.802 

1.000 

0.294 

0.924 

0.144 

0.746 

1.000 

0.193 

1.078 

0.081 

0.686 

1.000 

0.119 

1.232 

0.042 

0.624 

1.000 

0.068 

1.385 

0.020 

0.561 

1.000 

0.035 

1.539 

0.008 

0.498 

1.000 

0.016 

1.693 

0.003 

0.434 

1.000 

0.006 

1.847 

0.000 

0.373 

1.000 

0.001 

2.001 

0.000 

0.313 

1.000 

-0.001 

2.155 

0.000 

0.256 

1.000 

-0.001 

2.309 

0.000 

0.202 

1.000 

-0.001 

2.463 

0.000 

0.153 

1.000 

-0.001 

2.617 

0.000 

0.107 

1.000 

-0.000 

2.771 

0.000 

0.067 

1.000 

-0.000 

2.925 

0.000 

0.031 

1.000 

-0.000 

3.079 

0.000 

0.000 

1.000 

0.000 

ITAL  VERTICAL  MTFs 

cy/mr 

H_SYS 

H_PRE 

H_SPF 

0.000 

1.000 

1.000 

1.000 

0.200 

0.816 

0.967 

0.843 

4.712e+001  Hz 
4.714e+001  Hz 
3.101e4001  Hz 


W/(cm*cm*K) 

sqrt(Hz)/(cm*K) 
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0.400 

0.600 

0.801 

1.001 

1.201 

1.401 

1.601 

1.801 
2.002 

2.202 

2.402 

2.602 

2.802 
3.002 

3.203 

3.403 

3.603 

3.803 
4.003 


0.626 

0.450 

0.304 

0.193 

0.115 

0.064 

0.034 

0.016 

0.008 

0.003 

0.001 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


0.927 

0.881 

0.829 

0.772 

0.712 

0.649 

0.585 

0.520 

0.456 

0.394 

0.334 

0.277 

0.223 

0.174 

0.129 

0.090 

0.055 

0.025 

0.000 


0.674 

0.511 

0.367 

0.250 

0.162 

0.099 

0.057 

0.032 

0.016 

0.008 

0.004 

0.002 

0.001 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 


PREHLTER  VALUES  AT  NYQUTST 

horz  H_PRE(0.96)  =  0.73 1  vert  H_PRE(1 .25)  =  0.696 

SAMPLING  RATES 

horizontal  1.92  samples/mr 

vertical  2.50  samples/mr 

effective  2.19  samples/mr 

SENSOR  LIMITING  FREQUENCIES 


spatial  Nyquist 

horizontal  3.08  0.96 

vertical  4.00  1.25 

effective  3.51  1.10 

MRTD  3D  NOISE  CORRECTION  (AVERAGE) 

300  K  294  K 

horizontal  1 . 709  1 . 606 

vertical  1.709  1.606 
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MRTD  AT  300  K  BACKGROUND  TEMPERATURE 


cy/mr 

horz 

cy/mr 

vert 

cy/mr 

2D 

0.05 

0.154 

0.032 

0.05 

0.200 

0.042 

0.203 

0.042 

0.10 

0.308 

0.073 

0.10 

0.400 

0.098 

0.239 

0.050 

0.15 

0.462 

0.132 

0.15 

0.600 

0.183 

0.275 

0.060 

0.20 

0.616 

0.227 

0.20 

0.801 

0.325 

0.311 

0.071 

0.25 

0.770 

0.394 

0.25 

1.001 

0.582 

0.352 

0.084 

0.30 

0.924 

0.713 

0.30 

1.201 

1.064 

0.395 

0.099 

0.35 

1.078 

99.999 

0.35 

1.401 

99.999 

0.444 

0.117 

0.40 

1.232 

99.999 

0.40 

1.601 

99.999 

0.494 

0.139 

0.45 

1.385 

99.999 

0.45 

1.801 

99.999 

0.545 

0.165 

0.50 

1.539 

99.999 

0.50 

2.002 

99.999 

0.597 

0.195 

0.55 

1.693 

99.999 

0.55 

2.202 

99.999 

0.651 

0.231 

0.60 

1.847 

99.999 

0.60 

2.402 

99.999 

0.703 

0.274 

0.65 

2.001 

99.999 

0.65 

2.602 

99.999 

0.756 

0-325 

0.70 

2.155 

99.999 

0.70 

2.802 

99.999 

0.809 

0.385 

0.75 

2.309 

99.999 

0.75 

3.002 

99.999 

0.860 

0.456 

0.80 

2.463 

99.999 

0.80 

3.203 

99.999 

0.911 

0.540 

0.85 

2.617 

99.999 

0.85 

3.403 

99.999 

0.961 

0.640 

0.90 

2.771 

99.999 

0.90 

3.603 

99.999 

1.004 

0.758 

0.95 

2.925 

99.999 

0.95 

3.803 

99.999 

1.032 

0.898 

1.00 

3.079 

99.999 

1.00 

4.003 

99.999 

1.060 

1.064 

MRTD 

AT  294  K  BACKGROUND 

TEMPERATURE 

cy/mr 

horz 

cy/mr 

vert 

cy/mr 

2D 

0.05 

0.154 

0.033 

0.05 

0.200 

0.044 

0.203 

0.044 

0.10 

0.308 

0.075 

0.10 

0.400 

0.102 

0.239 

0.052 

0.15 

0.462 

0.137 

0.15 

0.600 

0.190 

0.275 

0.062 

0.20 

0.616 

0.236 

0.20 

0.801 

0.337 

0.311 

0.073 

0.25 

0.770 

0.409 

0.25 

1.001 

0.603 

0.352 

0.087 

0.30 

0.924 

0.739 

0.30 

1.201 

1.103 

0.395 

0.103 

0.35 

1.078 

99.999 

0.35 

1.401 

99.999 

0.444 

0.122 

0.40 

1.232 

99.999 

0.40 

1.601 

99.999 

0.494 

0.144 

0.45 

1.385 

99.999 

0.45 

1.801 

99.999 

0.545 

0.171 

0.50 

1.539 

99.999 

0.50 

2.002 

99.999 

0.597 

0.202 

0.55 

1.693 

99.999 

0.55 

2.202 

99.999 

0.651 

0.240 

0.60 

1.847 

99.999 

0.60 

2.402 

99.999 

0.703 

0.284 

0.65 

2.001 

99.999 

0.65 

2.602 

99.999 

0.756 

0.336 

0.70 

2.155 

99.999 

0.70 

2.802 

99.999 

0.809 

0-399 

0.75 

2.309 

99.999 

0.75 

3.002 

99.999 

0.860 

0.472 
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0.80 

2.463 

99.999 

0.80 

3.203 

99.999 

0.911 

0.559 

0.85 

2.617 

99.999 

0.85 

3.403 

99.999 

0.961 

0.663 

0.90 

2.771 

99.999 

0.90 

3.603 

99.999 

1.004 

0.785 

0.95 

2.925 

99.999 

0.95 

3.803 

99.999 

1.032 

0.931 

1.00 

3.079 

99.999 

1.00 

4.003 

99.999 

1.060 

1.103 

MDTD  AT  300  K  BACKGROUND  TEMPERATURE 
1/mr  MDTD 


0.20 

17.554 

60.173 

0.40 

8.777 

15.142 

0.60 

5.851 

6.803 

0.80 

4.388 

3.884 

1.00 

3.511 

2.532 

1.20 

2.926 

1.797 

1.40 

2.508 

1.354 

1.60 

2.194 

1.066 

1.80 

1.950 

0.868 

2.00 

1.755 

0.726 

2.20 

1.596 

0.620 

2.40 

1.463 

0.539 

2.60 

1.350 

0.476 

2.80 

1.254 

0.425 

3.00 

1.170 

0.384 

3.20 

1.097 

0.350 

3.40 

1.033 

0.321 

3.60 

0.975 

0.297 

3.80 

0.924 

0.276 

4.00 

0.878 

0.258 

4.20 

0.836 

0.242 

4.40 

0.798 

0.228 

4.60 

0.763 

0.215 

4.80 

0.731 

0.204 

5.00 

0.702 

0.194 
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MDTD  AT  294  K  BACKGROUND  TEMPERATURE 
1/mr  MDTD 


0.20 

17.554 

62.351 

0.40 

8.777 

15.690 

0.60 

5.851 

7.049 

0.80 

4.388 

4.024 

1.00 

3.511 

2.624 

1.20 

2.926 

1.863 

1.40 

2.508 

1.403 

1.60 

2.194 

1.104 

1.80 

1.950 

0.899 

2.00 

1.755 

0.752 

2.20 

1.596 

0.643 

2.40 

1.463 

0.559 

2.60 

1.350 

0.493 

2.80 

1.254 

0.441 

3.00 

1.170 

0.398 

3.20 

1.097 

0.363 

3.40 

1.033 

0.333 

3.60 

0.975 

0.308 

3.80 

0.924 

0.286 

4.00 

0.878 

0.267 

4.20 

0.836 

0.251 

4.40 

0.798 

0.236 

4.60 

0.763 

0.223 

4.80 

0.731 

0.211 

5.00 

0.702 

0.201 

FLIR92. . .  mitsh.l:  end  of  listing 
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